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Abstract. We consider a class of doubly nonlinear degenerate hyperbolic- 
parabolic equations with homogeneous Dirichlet boundary conditions, for which 
we first establish the existence and uniqueness of entropy solutions. We then 
turn to the construction and analysis of discrete duality finite volume schemes 
(in the spirit of Domelevo and Omnes [11]) for these problems in two and three 
spatial dimensions. We derive a series of discrete duality formulas and entropy 
dissipation inequalities for the schemes. We establish the existence of solutions 
to the discrete problems, and prove that sequences of approximate solutions 
generated by the discrete duality finite volume schemes converge strongly to 
the entropy solution of the continuous problem. The proof revolves around 
some basic a priori estimates, the discrete duality features, Minty-Browder 
type arguments, and "hyperbolic" L°° weak-* compactness arguments (i.e., 
propagation of compactness along the lines of Tartar, DiPerna, . . . ). Our 
results cover the case of non-Lipschitz nonlinearitics. 
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1. Introduction 

In this paper we consider degenerate hyperbolic-parabolic problems of the form 

{ dtu + div f{u) -diya{ V A{u)) = §, in Q := {0,T) x Q, 
u|t=o = uo, in Q 

u^O, on S = (0, T) X (917, 

where u : (i, x) e Q ^ R is the unknown function, T > is a fixed time, 17 C R*^ 
is a bounded domain with polygonal boundary dQ, and outward unit normal n. We 
consider the cases d = 2 and d = 3. The initial data wq : 17 ^ R are assumed to be 
a bounded measurable function, i.e., 

uoeL°°in), 

while the source § : Q — > R is assumed to be a measurable function for which 
S(i, •) e L°°{n) for a.e. t G (0,T) and /J" ||S(i, •)||L~(n) < oo; we abusively 
denote it by 

(2) §eL\0,T;L^iQ)). 
The function a : R^ R^ is taken under the form 

a(e) = 

where fc is a scalar function. The function a is assumed to be continuous and strictly 
monotone. We assume that there exist p G (1, +oo) and C > such that 

^ler' <fc(o <c|er-^ v^gr'^vio}. 
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In particular, the associated operator w ^ — div Vw|) Vw) is a Leray-Lions 

operator acting from Wg^'^(ri) to W~^'P (fl) with p' ~ A prototype example 

is the p-laplacian, which corresponds to k{^) — |^|^^^. 
We assume that the diffusion function A{-) satisfies 

A{-) is continuous and nondecreasing, normalized by A{0) = 0, 

while the convective flux function f(-) satisfies 

f = (fi, . . . , fd) : Q X R ^ R'' is continuous and normalized by f(0) = 0. 

We emphasize that the fluxes f, A are not necessarily locally Lipschitz continuous. 

Problems more general than (1), for which our results can be extended, will be 
discussed in Section 8. 

The class (I) of nonlinear partial differential equations includes several important 
particular cases. The hyperbolic conservation law 

dfU + div f (m) = 

is a special case of (1). The celebrated theory of L°° entropy solutions for scalar 
conservation laws in R"* was developed by Kruzhkov while the BV theory 
was set up by Vol'pcrt [77]. The extensions for the Dirichlet problem in bounded 
domains are due to Bardos, LeRoux, Nedelec [15] (for the BV setting) and Otto 
[09] (for the L°° setting). Note that the boundary condition is only verified in some 
generalized sense (see [15, 69, 65, 29, 7-3, 66, 79, 49, 67, 5]). 

Many other well-known partial differential equations (usually possessing more 
regular solutions) are also special cases of (1). Let us mention the heat and porous 
medium equations 

dtu = Au, dtu = Au™, m > 1, 
and more generally degenerate convection-diffusion equations of the type 
(3) dtu + d\Y^{u) ^ AA{u). 

Degenerate parabolic equations like (3) occur in theories of flow in porous media 
(see discussion and references [43]) and sedimentation-consolidation processes [27]. 

As other famous representatives of the class of equations that is considered 
herein, we mention the p-Laplace equation 

5tU = div(|Vu|^"^ Vw), p>l, 

which arises in the theory of non-Newtonian filtration. Also well known is the more 
general polytropic filtration equation 

dfU ~ div 



im — 1 

\u\ U 



\u\™' ^ u ) I , m,p > 1. 



A related class of equations consists of the so-called elliptic-parabolic equations 

dtb{v) = div a{v,\7v), 

where 5 : R — > R is continuous nondecreasing, and a(r, ^) : R x R^ R^ gives 
rise to a Leray-Lions operator. We refer to [4, 20, 70, 30, 6] and the references cited 
therein for more information on elliptic-parabolic equations. 

A chief goal of this paper is to propose and analyze a specific class of finite 
volume schemes for the problem (1). Note that finite volume schemes are well suited 
for approximation of equations in divergence form, such as (1). Discretization of 
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the aforementioned hyperbolic, porous medium, convection-diffusion, and elhptic- 
parabolic equations by finite volume methods is quite standard by now and often 
used in engeneering practice. We refer to [48, -il, -i, 44, 45, (if, 57, 08, 79, 49, 07, 
f 2, fO, f f , 42] and references therein for different convergence results and numerical 
experiments. For related works on linear elliptic problems, see [2, f , 57, 4f , 23, 58, 
50, 5f , 53, 52] and the discussion in Section 8. Alternative numerical approaches 
have also been investigated; here we only mention finite element schemes (see [30, 
fO] and references therein), kinetic schemes (see [f4, 22, 55] and references therein) 
and operator splitting schemes (see [43]). 

Having said that, we are not aware of any papers that construct convergent 
numerical schemes for mixed type equations of the generality considered herein. 
Indeed, they combine a number of difficulties such as nonlinear convection, doubly 
nonlinear diffusion, strong degeneracy, and shocks, which in turn necessitates the 
use of a suitable framework of discontinuous entropy solutions. Furthermore, in the 
absence of the Lipschitz continuity assumption on the convectivc flux f(-), the CFL 
condition docs not make sense; therefore wc have to discrctize the convectivc term 
with a time-implicit scheme. 

We begin by providing the entropy solution framework for (1); this is the topic 
of Section 2 and Appendix A. Due to the nonlinearity of f(-) and the possible 
degeneracy of A{-), the problem (I) will in general possess shock wave solutions, 
a feature that can reflect the physical phenomenon of breaking of waves. This 
is well known in the context of conservation laws. Also the boundary condition 
cannot be prescribed pointwisc on the whole boundary E when A is not strictly 
increasing. Due to this loss of regularity, it is necessary to work with weak solutions; 
moreover, to single out a physically relevant and unique weak solution, wc need to 
impose additional "entropy inequalities" , in the spirit of Kruzhkov [G3] . Early 
results on hyperbolic-parabolic equations were obtained by Volpert, Hudjaev [78]; 
see also [80, 82, 8f], [74], [28], [f9] and references cited therein, and [76], [37], 
[32]. entropy techniques for degenerate convection-diffusion equations like (3), 
which take into account both hyperbolic and parabolic features, were developed by 
Carrillo [29] for the homogeneous Dirichlet problem in bounded domains. Since 
then, many authors extended the Carrillo results in various directions (see e.g. 
[30, 59, GO, 73, 25, 40, 00, 02, 07, 49, 5, 13]). Some additional techniques are 
required for anisotropic diffusion problems, where a kinetic approach (see Chen, 
Perthame [35]) and an accurate cntropic approach (see Bcndahmane, Karlsen [17, 
18]) were developed in the few last years; see also Souganidis, Perthame [75] and 
and Chen, Karlsen [34]. In this paper, we use a variant of the Carrillo cntropic 
approach. Following the Tartar-DiPerna idea of measure- valued solutions and using 
the techniques of Eymard, Gallouet, Herbin [48], wc introduce a notion of entropy 
process solution for (1), and establish the related identification and uniqueness 
results. More exactly, we work with entropy double-process solutions arising in the 
particular context of discrete duality finite volume schemes. In Section 2, we show 
the existence result for (1) and state uniqueness; an adaptation of the standard 
uniqueness (and, more generally, contraction and comparison principle) proof is 
given in Appendix A. 

In Section 3 we construct discrete duality finite volume (DDFV) schemes for 
(1) in two and three spatial dimensions (some other schemes are briefly discussed 
in Section 8). We adapt the approximations used by Eymard. Gallouet, Herbin 
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[48] (see also [31, 79, 49, 67]) for the nonlinear convection term, and those used by 
Hermeline [57, 58], Domelevo, Omncs [ : ] and Andreianov, Boyer, Hubert [11] for 
the doubly nonlinear diffusion term. In 3D. we propose new DDFV schemes that 
possess convenient discrete duality properties. 

Our 3D scheme is a very particular case of the schemes introduced and studied 
numerically by Hermeline in [58]. In passing, we mention that different kinds of 
3D discrete duality schemes were constructed in [72, 3!)] and in [-SS]. Appendix B 
(see also [n, 7]) is devoted to an elementary reconstruction lemma which underlies 
our DDFV schemes in 3D. In contrast to [41, 11], we are led to penalize our DDFV 
schemes to ensure that the two approximations of A(u) actually converge to the 
same limit (see Section 3.4). The DDFV schemes constructed in Section 3 possess 
several convenient discrete calculus formulas that we collect in Section 4. Related 
consistency estimates and properties of the associated spaces of discrete functions 
are given in Section 5. The (few) available a priori estimates for the discrete so- 
lutions are collected in In Section 6. In the same section, the existence of discrete 
solutions is shown. Furthermore, we establish that, up to an error term in the equa- 
tion depending on the discretization parameter, discrete solutions can be considered 
as entropy solutions of (1). In Section 7 we prove that discrete solutions converge, 
as the discretization parameter tends to zero, to an entropy double-process solution 
that turns out to be the (unique) entropy solution of (1). It should emphasized 
that we obtain strong convergence of both convective and diffusive fluxes, in spite 
of the double nonlinearity of the problem (1). Section 8 contains references to some 
known finite volume schemes for nonlinear diffusion-convection equations, and dis- 
cusses the extension of our results to different generalizations of problem (1). 

2. Notions of solution and well-posedness 

As it was explained in the introduction, we need the notion of weak solution 
for (1) with additional "entropy" conditions. In order to use entropy conditions 
in the interior of Q and, moreover, take into account the homogeneous Dirichlet 
boundary condition on E, following Carrillo [29] we will work with the so-called 
"semi-Kruzhkov" entropy-entropy flux pairs (^7^,^^) for each c G R; they are 
defined as 

(z) = {z - c)+ , ?7^(z) = (z - c)-, 

q+(z) = sign+ (z - c) (f(z) - f(c)), q' (z) = sign~ (z - c) (f(z) - f(c)). 

By convention, we assign {r]^y{c) to be zero. Here (z — c)^ denote the nonnegative 
quantities satisfying z — c ~ {z — c)+ — (z — c)^; moreover, we use the notation 

sign^ (z — c) 
sign" (z — c) 

At certain points, we will also need smooth regularizations of the semi-Kruzhkov 
entropy-entropy flux pairs; it is sufhcient to consider regular "boundary" entropy 
pairs (?7^£,c|^g) (cf. Otto [09] and the book [()5]), which are W^'°° pairs with the 
same support as (?7^,qJ), converging pointwise to (?7^,q^) as e ^ 0. Specifically, 



(0, z<c, 

= (.c )'(.) = {«' 

1-1, z < c. 
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the functions 

signg^(z) — — niiniz''', e}, sign^(2:) = — max{— z^, — e} 

will be used to approximate sign^(-) ^ {vt )'(')■ 

In view of the monotonicity of A : R ^ R, the following definition is meaningful. 



Definition 2.1. For any locally bounded piecewise continuous function 9 : R 
we define (using, e.g., the Stieltjes integral) the function Ag : H —> R by 



(4) 



Aeiz) 



0{s)dA{s) 



The ensuing lemma shows that there exists a continuous function Ag such that 
Ag{z) = Ag{A{z)). We prove this lemma under rather strong assumptions, but 
they are still sufficient for our needs. 

Lemma 2.1. (i) Let 9,Aq be a couple of functions as introduced in Definition 2.1. 
Then there exists a continuous function Ag : A{R) R such that 

Ag{z) ^ Ag{A{z)), VzeR. 

Moreover, Ag is Lipschitz continuous. 

(ii) Assume additionally that 9 € W^'°°{R), and let {AP)p be a sequence of 
nondecreasing continuous surjective functions converging to A pointwise on R as 
p — + 0. Define Ag,Ag by (i) and (4) with A^ replacing A. Then Ag converges to 
Ag uniformly on compact subsets of A{R). 

Proof, (i) For b £ ^(R), we can define Ag by Ag{b) = Ag{z) for some z G A^^{b). 
If A{z) = A{z), then the measure dA(s) vanishes between z and z; thus 



Agiz)-Agiz) = / 9{s)dAis)=0, 

J z 

and Ag is well-defined. For all 6, 6 G ^(R), 

Ag{b) ~ Ag{b) ^ Ag{z) ~ Ag{z) = [ 9{s)dA{s), z £ A-\b) , z £ A-\b) . 

J z 

Consequently, 



Aeib)-Agib) 



< 



\Aiz)-A{z)\ = \\9\\, 



b-b 



(ii) Since the functions Ag arc monotone, by the Dini theorem it is sufficient to 
prove the pointwise convergence. By the same argument, the convergence of A^ to 
A is actually uniform on compact subsets of R. Take b £ A{R) and z £ A~^{b). 
Set b'' = AP{z); we have 6'' 6 as p — > 0. Using (i) and the integration-by-parts 
formula for the Stieltjes integral, we get 



A'g{b)^Ag{b) < A^gib) - A^gib") 

< \\9h^\b-bP\^ 



<2\\9\\^^\b-bf\ 



A^gib") ^ Ag{b) 
9{s) diAP{s) ~ Ais)) 



{AP{s) - A{s))9'{s)ds 



The right-hand side converges to zero as p — ^ 0. Thus the claim follows. 



□ 
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We have now came to the definition of an entropy solution. Here and in the 
sequel, denote {fc G R | ± fc > 0}, respectively. 

Definition 2.2 (entropy solution). An entropy solution of the initial-boundary 
value problem (1) is a measurable function u : Qt R satisfying 



(D.3) for all pairs (c, ?/;) G R^ x I?([0,T) x Vi), "0 > 0, and also for all pairs 
(c, v) e R X P([0, T) X n), i' > 0, 



For the convergence proof we need the notion of entropy double-process solutions; 
we adapt this notion from [48, 31, 54, 49], where entropy process solutions have 
been introduced for hyperbolic problems and degenerate parabolic problems with 
linear diffusion. This definition is based upon the so-called "nonlinear L°° weak-* 
convergence" property, which is well-known in the equivalent framework of measure- 
valued solutions developed earlier by Tartar and DiPerna: 

each sequence (up) bounded in L°°{Q) admits a subsequence such that 



where the function /x G L°°{Q x (0, 1)) is referred to as the "process function"; it 
is related to the distribution function of the Young measure. As usual, in (5) and 
elsewhere we do not bother to (re)label sequences. 

We remark that the reason for introducing in the definition below two different 
process functions fi,^J,*, both corresponding to the single unknown function u, is 
that it permits us to handle the double approximation of u by pairs u™ , in the 
framework of DDFV schemes (see Section 3). 

Definition 2.3 (entropy double-process solution). A triplet (^,/i*,w) of measur- 
able functions, with fi, ^J,* : Q x (0, 1) — > R and w : Q R, is called an entropy 
double-process solution of the initial-boundary value problem (1) if the following 
conditions are met: 

(DM) /i, ^* G i°°(Q X (0, 1)), w G LP(0, T; Wj'^(f^)), and 
A{^{t, x,a)) = w{t, x) = A(p* (i, x, a)), 
for a.e. (t, x,a) € Q x (0, 1). 



(D.l) u G L°°{Q) and w = A{u) G LP{0, T; Wl'^{Vl)); 
(D.2) for allipe V{[0,T) x n), 
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(D'.2) For all ^ e P([0, T) x Q), 

- / k{Vw)Vw Vi{;dxdt+ / Uoip{0,-)dx+ / §i(jdxdt = 0. 



(D'.3) For all pairs {c,ip) e R=t X V{[0,T) xn), ip > 0, and also for all pairs 
(c, v) e R X P([o, T) X rj), > 0, 



+ /Y ^((^,±)'(M) + (d-l)(77±)'(/x*))SVrfa;dtda>0. 



Remark 2.1. Since \7w — a.c. on {{t,x) G Q) | w{t,x) — A{c)} for any c G R, 
the term fc( Vu')VA^^±y(w) in the above definitions can be rewritten as 

(6) for any z G and also as sign''^(z« — A(c))a( Vw). 

The form used in (D.3) and (D'.3) is convenient for expressing the approximate 
entropy inequaUties at the discrete level; the equivalent form (6) is used in the 
uniqueness proof. Both forms arc exploited in the existence proof below. 

Remark 2.2. Let u be an entropy solution of (1). Then the triplet {jj., fj* ^w) 
defined by 

X, a) = p*(t, X, a) — u{t, x) for a.c. {t, x,a) G Q x (0, 1), 
w{t,x) = A{u{t,x)) for a.c. {t,x) G Q. 

is an entropy double-process solution of (1). 

Conversely, if {n,fj.*,w) is an entropy double-process solution of (1) for which 
fi{t,x,a) = ^*{t,x,a) = u(t,x) a.e. on Q x (0,1) for some function it : Q — > R, 
then this u is an entropy solution of (1). 

Note that in Definition 2.2, we have only considered a-independcnt data uo,/. 
In this case, the notion of entropy double-process solution is just a technical tool 
that permits to bypass the lack of strong compactness of sequences of approximate 
solutions. As a first illustration of this, wc pass to the limit in vanishing viscosity 
approximations (without BV estimates) to prove the existence of an entropy double- 
process solution such that fi = 

Theorem 2.1. Under the assumptions stated in Section 1, there exists an entropy 
double-process solution to the initial-boundary value problem (1) for which fi = fi* . 

Notice that the above result holds for any Lipschitz domain il in any space 
dimension. In passing, we also mention that the existence result of Theorem 2.1 
has recently been generalized by Ouaro and the authors [9] to the case of a triply 
nonlinear degenerate diffusion equation. 
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Proof. The proof is divided into several steps. 

(i) We approximate problem (1) by regular problems (1)^ with f, A replaced by 
fp, Ap such that fp, Ap, [Ap]~'^ are Lipschitz continuous on R and fp, Ap converge to 
f, A, respectively, uniformly on compacts sets as p ^ 0. 

Using classical techniques (cf. Alt, Luckhaus [4] and Lions [()4]), we can show 
that there exists a weak solution Up G LP(0, T; Wg^'^(ri)) to problem (1)^ in the 
following sense: 



(7) 



dtUp + divf p{up) =diva(VAp(wp)) + § 

in LP'{0,T;W-^^P'in)) + L\Q), Up\t^o ^ uo- 



Moreover, since fp o A~^ is Lipschitz continuous, the contraction property and 
comparison principle for weak solutions can be verified. It can be obtained either by 
the technique of Otto [70] (doubling the time variable) or using the theory of integral 
solutions and nonlinear semigroup methods, consult for example [30]. Besides, Up 
verifies the entropy formulation of Definition 2.2 with fluxes jp,Ap, where 77^ can 
be replaced by regular "boundary" entropies ry^^, whenever we prefer to do so. 
(ii) We claim that the following quantities are uniformly bounded in p: 

• l|Mp||^oo(n) and ||^p(wp)|lip(o,T;Wo'P(n))' 

• space translates of Ap{up) in L^{Q) (consequence of previous estimate); 

• time translates of Ap{up) in L^{Q). 

Indeed, for the first point consider the function 



M{t) ||uo|Il=c,(s:2) + ||§(t, ^[^-(o) dr, 

which is a solution of (l)p with .x-constant data ||uo||^oo(q), II§(^i ■)llL=°(n)- "^^^ 
comparison principle mentioned in (i) ensures that a.e. on Q, 

-~M{T) < -M{t) < Up{t,x) < M{t) < M{T). 

Next, we employ j4p(up) as a test function in (7). The product between dtUp and 
Ap{up) is handled using the usual chain rule argument (see, e.g., [4, 70, 30]), where 
the relevant duality is between the space E := LP(0, T; Wp^'^(17)) n i°°((3) and 
the space Lp'(0,T; W'^'P' {Q)) + L^{Q) C E*. Here we are also exploiting the 
bound on fp(up) in a straightforward fashion to treat the term fp(wp) • VAp{up); but 
notice that using the Green-Gauss trick (16) below, we can supply a finer analysis 
of this term. 

For the third bullet point, we first use (7) to get, for a.e. t, i + A G (0, T), 

i-t+A P 

{up{t + A) - up{t)) e = / / [ (-fp(iip) + a( VApiup))) • Ve + § e ] 

foraUC e W^'''{n)nL°°{n). Taking^ = Ap{up{t + a)) - Ap{up{t)) and integrating 
in t, using the two previously obtained estimates, we deduce that 



(8) // |up(t + A)-Up(t)||^p(up(t + A))-Ap(up(t))| <Const|A| 



Q 



Now, let TT be a (common for all p) concave modulus of continuity for Ap on 
[—M{T), M{T)], n be its inverse, and set Til{r) ~ rn(r). Let tt be the inverse 
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of n. Note that tt is concave, continuous, and tt (0) = 0. Set v{t, x) = Up(t + A, x) 
and y(i, x) = Up{t, x). We have 

Ap{v) - Ap{y)\ = l^n(n{\Ap{v) - Ap{y)\) 

<\Q\n(^^J^Il{\A,{v)-Ap{y)\ 

Since \Ap{v) - Ap{y)\ < tt{\v - y\), we have U{\Ap{v) - Ap{y)\) < \v - y\ and 
II(|^p(z;) - Ap{y)\) = U{\Ap{v) - Ap{y)\)\Ap{v) - Ap{y)\ 
< \v-y\\Apiv)-Apiy)\. 

Therefore, estimate (8) impHes 

\Ap{up{t+A, x)) - Ap{up{t, x))\ 

Q 

(9) <\Q\n(^jLJ^\y^y\\Ap{v)--Ap{y) 

1 



l^|J(A)j <CniCA) =:uoa{a), 

where uja G C(R+,R+), wa(0) = 0. 

(iii) Thanks to the estimates in (ii) and standard compactness results, there 
exists a (not labelled) sequence p ^ such that 

• Wp = Ap{up) converges strongly in L^(Q) and pointwise a.e. on Q; 

• Vwp converges weakly in LP{Q); 
a{ Vwp) converges weakly in (Q) to some limit x; 

• Up converges to : Q x (0, 1) G R in the sense of (5). 

Let us introduce the function 

(10) u{t,x)~ / ii{t,x,a)da, for a.e. {t,x)&Q. 







Thanks to the convergence of Ap to A, we can identify the limit of Wp(-, •) with 
A{^{-, •, a)) da. Moreover, since Wp is converging strongly, A(/x(-, •, a)) is actually 

independent of a G (0, 1) and equals A{u{-, •)). Using distributional derivatives, we 

also identify the limit of Vwp with VA{u). 

(iv) We have now come to the main step of the proof, namely to improve the 

weak convergence of \7Ap{up) to strong convergence, and to identify the weak 

limit of a( VAp{up)) with a( VA{u)), where u is defined in (10); of course, the chief 

difficulty comes from the lack of strong convergence of Up. 

We begin by specifying the test function in (7) as Wp C, yielding 

(11) / {dtup,wp0~ [ fp{up)- \7wpC+ f a{\/wp)- Vwp C - / ^WpC^O, 



Q JQ 



where Wp = Ap{up) and C G T>{[0,T)) is nonincreasing with C(0) = 1. Next, we 
pass to the limit into the weak formulation (7), obtaining 

j dtu + div f(/^) da = div x + § 

I in LP'iO,T;W-^^P'iQ)) + L\Q), «p|t=o = "o- 



(12) 



FV FOR DOUBLY NONLINEAR DEGENERATE EQUATIONS 



11 



In (12), wc take w C sls test function, where w = A{u), u is defined in (10), and C, 
is as specified above. The resuh is 

(13) / {dtu,w0- [ [ m- ^ujC+ [ X- ^A{u)C- [ §wC = 0. 
Jo ^ JQJa JQ JQ 

In order to later use the Minty-Browder trick, we shall combine (13) and the 
"p 0" limit of (11) to conclude the validity of the following inequality: 



(14) / X • VA(w) > liminf / a{\7wp) ■ Vwp. 

Jq p^'^ Jq 

A crucial role is played by the following calculation, which reveals that the 
lack of strong convergence of fp{up) is not an obstacle. Indeed, a componentwise 
application of Lemma 2.1 (i) yields the existence of a Lipschitz continuous vector- 
valued function such that 



(15) / K-^)dAis)=MA{z)). 

Jo 

Hence, by the chain rule and the Green-Gauss formula, we can calculate as follows: 

/ /'f(M) • vA{u) = / ^\{^,) ■ VA{^,) = f' r f dwMA{^,)) 

JqJo jqjo Jo jo Jn 

I A^{A{u))-n^Q, 
J on 

because for a.e. a G (0, 1), 

A{p{; ; a)) = A{u{; •)) e LP{0, T- <'^(f^))- 
By similar (simpler) arguments and Up e i^(0, T; Wq"^{Vl)), we also have 

(16) / fp(up) • VAp{up) = / / div f / ' ]p{s)dAp{s)) = 0. 
JQ JO Jn \Jo / 

Consequently, we can make I2 and /2,p (for each p > 0) vanish. 

Next, let us prove that Ii < limp^o-^i.p- As above, the duality products 
{dtUp,Ap{up)), {dtU,A{u)) are treated via the chain rule argument (cf. [4]). Set 
B{z) = /g A{s) ds, Bp{z) = /p Ap{s) ds, and note that these functions are convex. 
Also, Bp ^ B uniformly on compact subsets of R. With the help of Jensen's 
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inequality, 

{dtu,A{u)0 = - f B{u)C- f Biuo) 
/o JQ Jn 

< f f Bifiit, X, a)) da i^O- I B{uo) 
JQ JO Jn 



lim f- / SpK)C'- / B, 
V Jq Jn 



("o) 



T 



= lim / {dtUp,Ap{up)(} = lim h^p. 

P^O P^O 

Finally, it is clear that I^^p — *■ /s as p — > 0. Letting ( tend to 11[o,t)i the desired 
inequality (14) follows from subtracting the "p — > 0" limit of (11) from (13) and 
the above calculations. 

Starting off from (14). we can use the Minty-Browder trick (see, for example, 
[64, 24, 4, 21] and the proof of Theorem 7.1 in Section 7) to deduce that 

(17) a( Vwp) - a{ VA{u)) weakly in L^' (Q) as p ^ 0. 

Thus X ~ a{\7A{u)). Simultaneously, from the strict monotonicity of a(-) we 
deduce that, firstly, the convergence in (17) also takes place a.c. in Q; secondly, 
that (14) actually holds with an equality sign. Next, we consider the functions 
5p := a( Vwp) • Vwp > and g := a( \7A{u)) ■ VA{u) > 0, and observe that 



gp ^ g a.e. in Q, 9p ^ 9 as p ^ 0. 

JQ JQ 

Hence, we deduce that a subsequence of {gp)p converges to g strongly in L^{Q), 
cf. [24], [21, Lemma 5], [42, Lemma 8.4]. Due to the coercivity of a(-), (| Vwp\P) ^ 
is equi-integrable, so the Vitali theorem yields the strong convergence of ViUp, 
along a subsequence if necessary, to a limit already identified as Vw, w = A{u). 

(v) By (12), we readily conclude that (^, fi, w) verifies (D'.2). Now we can pass 
to the limit in the entropy inequalities corresponding to (l)p and deduce (D'.3). 

Let us first show that VA'' j- s,(u'p) converges weakly to VA, ± y{w) in LP{Q). 
By Lemma 2.1 (i). A'' ± (•) are uniformly Lipshitz continuous functions. Thus 

\Vc,e) 

VA'' ± y(''^p) are uniformly bounded and weakly compact in LP{Q). Moreover, 
A'' ± v(^p) converges to A, ± y(i/;) by Lemma 2.1 (ii) and because of the pointwise 

\Vc,e) ^ lc,e J 

convergence of Wp to w. Using the distributional convergence, we eventually work 
out our claim. 

Now note that Hp > 2, then k is continuous. By the last result of (iv), we 
can assume without loss of generality that k{ \7wp) converges to fc(Vu') a.e. in Q. 
Moreover, (fc(Vii;p)) is bounded in Lp^{Q), since (Vu>p) is bounded in LP{Q). 
Applying the Egorov theorem and Holder's inequality with exponents p',p in the 
product (fc( Vwp) Vt/i) ■ VAj-^^i y{wp), we deduce that 

(18) lim / k{Vwp)yA,± y{wp) • W / fc( Vw) VX + y{w) ■ Vip. 
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li p < 2, we fix a small S > and truncate fc(-) in the (S-neigbourhood of the origin 
(if k{-) is replaced by ks{-) = min{A:(-), min|^|<5 the argument used for p > 2 
applies), and we analyze separately the set | | Vw''{t,x)\ < S}. On this set, 

kiVwp) I VI(„±^),K)| <C\\i7^tJ\Uywpr' < Const -^^--i, 

uniformly in p. To conclude that (18) still holds, we first pass to the limit as p ^ 
for a fixed 5 > 0, and then send J — > 0. 

Let us take regular "boundary" entropy pairs r]^^ such that [rj^^)' approximate 
{rf^y (extended by zero at the point c, by our convention), pointwise a.e. in R as 
£ ^ 0. We use (18) to pass to the limit in the entropy inequality corresponding 
to (l)p. We pass to the limit in the remaining terms in this entropy inequality 
using the continuity of t?^^, q^^, (?7^e)' and the nonlinear L°° weak-* convergence 
property (5). Finally, we pass to the hmit as £ — > 0, rewriting VA^^± y(w) as 
ivtsYi'^) (consult Remark 2.1) and using the Lebesgue dominated convergence 
theorem and the pointwise convergences of '7^e, q^^, (jy^^)'. The passage to the 
limit in the weak formulation is similar. 

(vi) Wc conclude that (/i, fi*,A{u)) is an entropy double-process solution of (1) 
such that 1^1* ~ fi. □ 

Given Theorem 2.1, the uniqueness of an entropy double- process solution can be 
established using Kruzhkov's method, along the lines of Carrillo [2!)]. 

Theorem 2.2. Suppose the assumptions stated in Section 1 hold. Let (^,/i*, w) he 
an entropy double-process solution of the initial-boundary value problem (1). Then 
it is unique. Moreover, there exists a function u G L°°(Q) such that 

fi{t, X, cx) = u{t, x) = /i*(i, X, a) for a.e. (t, x,a) ^ Q x (0, 1). 

We refer to Appendix A for a sketch of the proof. 

Theorems 2.1 and 2.2 as well as the arguments of Appendix A imply 

Corollary 2.1 (well-posedness). Under the assumptions stated in Section \, there 
exists a unique entropy solution of the initial-boundary value problem (1). Let u 
and V be two entropy solutions of (1) with initial data u\t=o — uq & L°°{Q) and 
v\t=o — vq £ L°°{Q) and source terms § and T of the kind (2), respectively. For 
a.e. t € (0,r), we have 

I {u{t,x) - v{t,x))^ dx < [ (uq-vq)^ dx+ [ [ (§~T)+. 
Jn Jn Jo Jn 

Consequently, if uq < vq a.e. in and § < T a.e. on Q, then u < v a.e. in Q. 

Finally, if uq = vq a.e. in Q and § ~ 7 a.e. on Q, then u — v a.e. in Q. 

The upcoming sections are concerned with the construction of finite volume 
schemes for which the corresponding discrete solutions converge to the unique en- 
tropy solution of (1) as the discretization parameter (mesh size) tends to zero. The 
convergence proof will attempt to mimic the proof of Theorem 2.1. 

3. Discrete duality finite volume (DDFV) schemes 

Let n be a polygonal (respectively, polyhedral) open bounded subset of R'', 
d = 2 (respectively, d = 3). In what follows, we introduce most of the notation 
related to DDFV schemes; each piece of new notation is given in italic script. 
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3.1. Construction of "double" conformal meshes. 

• A partition of is a finite set of disjoint open polygonal (respectively, poly- 
hedral) subsets of such that is contained in their union, up to a set of zero 
d-dimensional measure. 

Following Hermeline [ ], Domelevo, Onincs [4 i] and Andreianov, Boyer, Hubert 
[11], we consider a DDFV mesh which is a triple T = (OJl, Wl* , &) described below. 

• We let 271 be a partition of into triangles (respectively, tetrahedra); a more 
general case is discussed in Section 8^. We assume that the mesh satisfies the 
Delaunay condition (see, e.g., [4ti]); for simplicity of the representation, the reader 
may assume that each triangle (respectively, tetrahedron) contains the centre if its 
circumscribed circle (respectively, ball). We assume in addition 

, , if d = 3, each face of each tetrahedron of StJl 

contains the centre of its circumscribed circle. 

Although the definition of the scheme does not require condition (19) (see Re- 
mark 3.1 below), we do need this condition in order to deduce the discrete entropy 
inequalities and to prove that the scheme converges. 

Each control volume K £ fJJl is supplied with a centre Xk that we choose to be 
the centre of the circle (respectively, ball) circumscribed around k. We call 
the set of all edges (respectively, faces) of control volumes that are included in 
90. These edges (respectively, faces) are considered as boundary control volumes] 
for K € dVn, we choose the middle of k (respectively, the centre of the circle 
circumscribed around A') for the centre Xk- We denote by Wl the union miUdWl. 
We call vertex (of 931) any vertex of any control volume k E 9K. 

• (see Figure 1) We take 93T* as the partition of Q into dual control volumes if*, 
supplied with dual centres Xk* , such that vertex of 9JI and k' is the subset 
of points of n that are closer^ to Xk" than to any other vertex of 93T. In other 
words, 9Jl* is the Voronoi mesh constructed from the vertices of 231. If Xk* G 
we say that k' is a dual control volume and write k* S 93T*; and if Xk* G 

we say that k* is a boundary dual control volume and write k' G Thus 
SOT* = 9Jt* U dVn* . We call dual vertex (of SPT* ) any vertex of any dual control 
volume K* g Q3T*. Note that by the choice of x^, the set of centres coincides with 
the set of dual vertices, and the set of vertices coincides with the set of dual centres. 
In other words, S!3T and S!3T* are finite volume meshes that are dual each one to the 
other. 

• We call neighbours of K, all control volumes l G SUT such that k and l have a 
common edge (respectively, common face). The set of all neighbours of k is denoted 
by J^{k). Note that if l S J^{k), then k G in this case we simply say that k 
and L are (a couple of) neighbours. 

• (see Figures 1 and 2(b)) If k and l are neighbours, we denote by the 
interface dx n Ol between k and l. The set of all interfaces is denoted by £. 



In particular, in the two dimensional case we can partition f2 into polygons that admit a 
circumscribed circle. In the three dimensional case , we can partition H in polyhcdra that have 
triangular faces and admit a circumscribed ball. 

"^in order to avoid pathological situations which could appear in non-convex domains, e.g., in 
domains with cracks, here the distance between two points x,y ol is understood as the length 
of the shortest path which connects x with y and which lies within f2 
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Figure 1. 2D primal and dual meshes; diamond (subdiamond) . 



• In the same way, we denote by N'*{k*) the set of (dual) neighbours of a dual 
control volume k' , and by if*|L*, the (dual) interface Ok* HOl* between dual neigh- 
bours K* and i*. The set of all dual interfaces is denoted by £* . 

• (see Figure 2) The meshes Wl and induce partitions of 51 into diamonds 
and subdiamonds. Let us describe them separetely for d = 2 and d = 3. 

For d = 2 (see Figure 2(a)), if if , l e 9Jl are neighbours, then there exists a 
unique couple of dual neighbours {k', l*} such that the interface k|l is the segment 
with summits Xk* and Xl* ■ Then the quadrilateral £>x*i,* which is either the union 
(if Xk,Xl lie on different sides from _k|l) or the difference (if Xk, x^ lie on the same 
side from i<jL) of the triangles XkXk'X[^», x^Xk'Xl* is called a diamond; it is also 
unambiguously denoted by d^'^. 

For d = 2, every diamond is also called a subdiamond; the subdiamond which 
coincides with a diamond d'^^'^ is denoted by s^l^,. 

For d = 3 (sec Figure 2(b)), ii k,l G Wl arc neighbours, then there exists a 
unique triple of dual neighbours {k' , l' , m'} (which arc neighbours pairwisc) such 
that the interface k\l is the triangle with summits Xk*, x^* and Xm*- Then the 
polyhedron D^l'^.j^j. which is either the union (if Xk,Xl lie on different sides from 
i^) or the difference (if Xk,Xi^ lie on the same side from k\l) of the pyramids 
X/cX/^^Xi.Xjj*, X i^x K* X X is called a diamond; it is also unambiguously denoted 
by D^'^. Each diamond is split into three subdiamonds; e.g., the subdiamond sjjj.^* 
is the convex hull of Xk, Xk- , a^i,, Xi^* . 

We denote by 3?,© the sets of all diamonds and the set of all subdiamonds, 
respectively. Generic elements of S),(5 are denoted by respectively. 

Remark 3.1. If we drop condition (f9), the orthogonal projection of Xk (which 
coincides with the projection of .Xi) on may not be contained within A'jL. To 
cope with this problem, one could consider subdiamonds of signed volume, not 
necessarily contained within the corresponding diamonds. Up to a permutation 
of the subscripts a-*, i,*, m*, we have instead of the decomposition d^,*^*^^. = 
^M^K*-! ^^^^ decomposition d^,j^,^j, — \Sk',l* ^l',m'J \ ^m*k* 'j 
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in this case the volume of s^^.^^, will be taken with the sign "minus". Under this 
convention, Lemma 3.1 below holds true, so that formulas (22), (23)-(25) below still 
yield consistent discrete gradient and discrete divergence operators which enjoy the 
discrete duality property [7]. But the discrete entropy dissipation inequalities of 
Proposition 4.2 would fail, which undermines the subsequent convergence analysis. 

• For all bounded set E C R**, set diam [E) = sup^. — x|. 

• We denote by the measure of an object E in its natural dimension (i.e., 
the c?-dimensional measure, if is a control volume, a dual control volume, a 
subdiamond or a diamond; and the (d— l)-dimensional measure, if E is an interface 
or a part of an interface). According to Remark 3.1, for the definition of the scheme 
we could drop (19), in which case for a subdiamond S^i^^. such that S^i^^. Dd^'^ = 
its volume is taken with the sign "minus" . 

3.2. Mesh parameters and regularity of meshes. 

• We define the size of the mesh by size(X) = '^^'^Eemuloi^u^ diam [E). 

• Following [11], we call the maximum among 



1/ */ ^^ (diam(K))'' (diam (k*))'* 

max card(A^ (^*)), max , max- 



rriK ^ "rriK* 



max 

Kn 



/diam (k) diam (d) \ /diam(K*) diam (r>) 

D^(Z> \ diam [d) diam [k) / k* no^id \ diam (z>) diam {k*' 



(where the maximums are taken over all k G SUT, k* G SOI*, r> G 3?) the regularity 
constant of the mesh and we denote it by rcg(X). Roughly speaking, this constant 
controls the ratio of dimensions of neighbouring control volumes, diamonds and 
dual control volumes, as well as the proportions of each volume. 

In all the discrete estimates and convergence results stated below, we require 
the family of meshes {'Zhjh to have regularity constants reg(Xft) that are uniformly 
bounded in h. In the sequel, whenever there is a dependency of various constants 
on reg(X), we tacitly assume that this dependency is increasing. 




(a) 2D (sub)diamond. (b) 3D primal volumes, diamond, subdiamond. 



Figure 2. Diamonds and subdiamonds. 
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3.3. Discrete gradient and divergence operators. Diamonds permit to define 
the discrete gradient operator, while subdiamonds permit to define the discrete 
divergence operator (see (20), (21), (22) and (23), (24), (25) below, respectively). 
Both arc needed to discrctize the second order "diffusion" operator in equation (1). 
But first we need to introduce some more notation. 

• (see Figure 3) For a subdiamond s = s'^l\,, wc denote hy a = as, <J* = <jI 
the (parts of the) interfaces s D and s fl k'\l', respectively, and by i^s, i^t, unit 
normal vectors to CTs and cr* , respectively (their orientation is chosen arbitrarily) . 

• For a diamond d = d^'^, we denote by Proj^, Proj* the operators of orthogo- 
nal projection of R'' on the subspaces < x^Xj^ > and on < x^Xj^ >"*-, respectively. 
One should note that wc have < i^s >=< XkX^ > and < ><Z< x^Xj^ >^ for all 
s e & such that s C d. 



• For a couple of neighbours if , l £ Wl, denote by dxt, dK.K\L, and v^^l the 
distance between Xk and , the distance from Xk to i<|L, and the unit normal 
vector to i<|L pointing from if to l, respectively. More generally, if if G 9Jl, then 

denotes the exterior unit normal vector to Ok. In the same way, for neighbours 
if*,L* e Wfl* we define dK*L', dK'.K*\L', and i^k'l'', for A'* G we define Vk*- 

Remark 3.2. Note that by construction both meshes 371, WT are conformal (or- 
thogonal) in the sense if [4^]); combined with the Delaunay condition, this means 
that ly^.L ■ XkXz, — dxL, i'k'l* ■ Xk*Xl* = dK*L* for all neighbours k,l and k'^l*, 
respectively. 

The conformity property is particularly important for our framework imposed 
by the possible degeneracy of the diffusion term and the presence of the hyperbolic 
convective term. On the other hand, if this term is dropped, non-conformal double 
meshes can be considered for c? = 2 (see [57, 41, 11]) and d = 3 (see [72, 39, 58, 8, 
7, 38]) within the variational framework. 

• A discrete function on f2 is a set w"^ — (w™, u™') consisting of two sets of real 
values w™ = (u'k)^^™ and w™* = {wK*)K*G9n' ■ The set of all such functions is 
denoted by R'^. 

A discrete function on is a set 



consisting of four sets of real values 



interface 




Figure 3. Notation in a subdiamond (2D and 3D). 
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The set of ah such functions is denoted by R"^. In case all the components of w^™ 
and of ui^™' are zero, we write w'^ G RJ . 

• A discrete field on $7 is a set JF'^ = (To) of vectors of R''. The set of all 
such functions is denoted by (R'')'°. 

• On the set R'^ of discrete functions w'^ on fi, we define the discrete gradient 
operator V^[-] by 

(20) : e R^^ V^w^= (V^w^)^^^ G (R'')^ 

where V^w^ is the discrete field on fl with values 
for d^2: 

^ ^ Wi - Wk Wr,* - Wk' , K L K.L 

(21) Vou;^ = — v^^r.^ ^ i^/c*!.* for u = _d'^-^ = s^.,,.; 

for d = 3: 

VoW = 1'k,L H m„K.L l^K^L* 

dKL rrio \ ^K'^L' dK'L' 

Wm' - w,, w^. - w„ 



(22) 



+ mgK,L iyj:^*M*+nigK.L J^M*Jf' 

K.L ^K.L I I ^K.L I I „K.L 



for D ^ D^i^,,,. = sit,, U sU":,, U 5 



Remark 3.3. Formulas (21) and (22) have the following common meaning. The 
vector Vdw'^ is the unique element of R'^ such that Proj ( VDif"^) = ^k'l' ■ 

^K*L* ' 

Further, for d ~ 2, Proj* ( Vjjw'^) is the gradient of the (unique) affine function 
on the interface -k|l (which is a segment with summits Xk*,Xl*) that takes the 
values Wk'iWj^* at the points x^^* and x^*, respectively. Similarly, for d = 3, 
Proj* ( VdW'^) is the gradient of the (unique) affine function on the interface k\l 
(which is a. trictngle with, sinnmits , x^^* , Xj^^.j* ) that takes the values ,u'm* 
at the points Xk* ,Xl* ,Xm* , respectively. 

Thus, the primal mesh fDT serves to reconstruct one component of the gradient, 
which is the one in the direction XkxI. The dual mesh 2Jl* serves to reconstruct 
the {d—\) other components which are the components in the (d— l)-dimensional 
hyperplane containing and is orthogonal to x^xl. 

The first and second assertions of Remark 3.3 are evident. Note that formula 
(21) easily generalizes to quite arbitrary non conformal double meshes (see [11, 
Lemma 2.4]). The third assertion is a direct consequence of the 2D reconstruction 
result of Lemma 9.6 given and proved in Appendix B (see also [iS, 7]). 

Remark 3.4. The discrete gradient is exact on affine functions. More precisely, let 
n be a diamond (n = d'^',^^, , \i d ^ 2; d ^ d'^:,^^, , if d = 3). Let w{x) := WQ+r-x, 
Wo, r G R'', be an affine function. If w'^ is a discrete function with values 

Wk = w{xk), Wl = w{xl); 

Wk* = w{xk*), Wl* — w{x^* ) (and Wm* — w{x^,* ) if c? = 3), 

then V ow'^ ~ r = Vw. This property follows by a straightforward comparison of 
the formulas (21) and (22) for the discrete gradient with the reconstruction formulas 
of the next lemma. 

Lemma 3.1. Consider d = d'^'^ G 2). With the notation above, for all r G R'' 
one has the following reconstruction properties: 
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ford = 2, r = [r ■ Vk,l) Vk.l + {r ■ Vk^i^.) Vk',l'] 
for d = 3, 

Proof. For d = 2, the claim is a straightforward consequence of the conformity 
of the meshes (see Remark 3.2); Vk,l, J^k'l* form an orthonormal basis of R^. 
When d = 3, the claim follows from the orthogonality of v^ ^ to _KjL and from the 
2D reconstruction property of Lemma 9.6 (cf. Appendix B) applied in the plane 
containing _k|l- D 

Remark 3.5. The fourth assertion of Remark 3.3 indicates possible generalizations 
to the multi-dimensional case. Unfortunately, it can be shown that if d > 4, the 
direct generalization of the reconstitution formula of Lemma 3.1 holds only for 
meshes 2Jl with very special geometries, such as the uniform simplicial meshes 
(see Remark 9.2, which has to be combined with an induction argument on the 
dimension d in order to link the weighted projections on the edges appearing in 
Lemma 3.1 with the weighted projections on the faces appearing in Lemma 9.6). 

• For s € & such that s C d with d G 2), we assign Vgw'^ = V^u'^. More 
generally, if T'^ is a discrete field on fi, we assign J^s = for s C d. 

For K G 9Jl, we denote by v{k) the set of all subdiamonds s G © such that 
Kds ^ 9- In the same way, for k* G SDT* we define the set v*{k') of subdiamonds 
intersecting k'. 

• On the set (R'')® of discrete fields J-"^ , we define the discrete divergence oper- 
ator div "^[-J by 

(23) div^ : G (R'')^ ^v^ ^ diY^[T^] G R^, 
where the discrete function = ("f ™, on ft is given by 

(24) w™ = {vkJkgwi with Vk = — , m^sJ^s ■ where = Vk\s] 

(25) W™* = (Wk*)a-*6™*i Wk*=^— , ma'Ts-VK', Vk'^VkAs- 

In (24), (25) for s given, — i^k\s denotes the restriction on Us of the unit normal 
vector Vk to dK exterior to if; therefore it means the one of the vectors i/g, —Vs 
that is exterior to k (see Figure 3). Similarly, v^- ~ v^- \ s is the one of the vectors 
V*, — 1/| that is exterior to k* . 

In fact, formulas (24), (25) can be conveniently expressed in terms of vector 
products involving the discrete field Ts and specific geometric objects depicted in 
Figure 3 (see [S, 7]). 

3.4. Penalization operator. On the set R'^ of discrete functions w'^ on fJ, we 
define the operator of double mesh penalization by 



-.w^ e R'^ = V'^[w'^] G R^, 
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where the discrete function v'^ = (w™,^™*) on il is given by 
(26) = (?;k)kG3« with ^ (d-l) . ^ — m^nx* (w^ - w^- ); 



sizc(T) m,^ 

Got* 



(27) v""' ^ {vK*)K'e^' with Vk' ^ -. — 7^ y^m,tnK'{wK' - w^). 

size(T) TO^. 

The penahzation is needed in order to ensure (without using the strong convergence 
of V^ui'^, cf. the proof of [11, Theorem 5.1]), that the two components of a discrete 
"double" function w'^ converge to the same limit. 

Remark 3.6. The choice of penalization operator we propose here is just the 
simplest possibility. In (26), (27), the difference {wk* - Wk) could be replaced by 
-Wk), which seems more natural with respect to the assump- 
tions on o; the power of size(T) in the denominator can be chosen arbitrarily. The 
convergence of the scheme would remain true. The question of optimal choice of 
the penalization operator is beyond the scope of this paper. 

3.5. Discrete convection operator. Let f : R ^ R*^ be continuous. Denote by 
LUi\i{-) a modulus of continuity of f on [— il/, AI], i.e., a continuous concave function 
on [0, M] with ujm{0) = and 

, , llf(«)-fWII <^A/M- 

a,bl£[~M,AI],\a~b\<r 

Note that we can always choose com strictly increasing, upon replacing lo]\j by 
COM + Id if needed. 

Following Eymard, Gallouct. Hcrbin [48], we now define discrete convection 
fluxes, separately for each of the meshes Wt, OJl*. This will allow to discretize 
the convectivc part of equation (1). 

• Let k\l G £. To approximate f{u) ■ Vk.l by means of the two values Uk, that 
are available in the neighbourhood of the interface k\l, let us use some function Qk^^ 
of the couple {uk, Ul) G R^. More exactly, take a collection of numerical convection 
flux functions (gK,L)K\L£e, 9k,l G C(R^,R); with the following properties: 

' (a) gK,L{',b) is nondecreasing for all & G R, 

and gK,L{ci, •) is nonincrcasing for all a G R; 
(b) gKL{a, a) = f(a) • Vk^l for aU a G R; 

(28) \ (c) gK,L{a, b) = —gK,L{b, a) Va, 6 G R, for all neighbours k, l G SUT; 
(d) gK,L has the same modulus of continuity as f, i.e., 

there exists C independent of k]l such that Va, b,c,d E [— M, M], 
\gK,L{a,b) - gK,L{c,d)\ < C {u!M{\a - c\) + ujM{\b - d\) ). 

These assumptions (see [48]) are by now standard. Note that the assumption 

(28) (d) usually states that f,gK,L are Lipschitz continuous, with the same Lipschitz 
constant; here, we adapt it to the case of general continuous function f. 

Note that (28) (b) and (c) are compatible. Also note that the consistency re- 
quirement (28) (b) together with the Green-Gauss formula imply 

(29) rnKiLgKhia, a) ~ Ha) I i'k = for all a G R, for all k G 9Jl. 
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Practical examples of numerical convective flux functions can be found in [4n]. 
These include the Godunov, Lax-Friedrichs, Engquist-Osher and Rusanov fluxes as 
particular cases. 

• Numerical convective flux functions gn'L*, k'^' G E* , are defined similarly. 

• On the set R'^ of discrete functions on fJ, we define the operator (divcf)'^[-] 
of discrete convection by 

(divcf)^ -.v^ elC ^ = (divcf)'^[u^] e R^, 

where the discrete function v'^ = on f2 is given by 

= (wrOa'g™ with u^, = — mKiL9K.L{uK,UL); 

3.6. Projection operators and test functions. 

• On L^iVl), we define the mesh projection operator on the space of discrete 
functions on D. by 

: S e L^{Vl) ^ = P^[S] e R^, 
where the discrete function S'^ = (S™, S™') on Vl is given by 

S™ = (§R-)A-e™ with Sk = — / 2>{x)dx\ 



(30) 

§™* = (S7,0a-g™- with = / S(a;)da;. 



rriK' Jk 

• For a sufficiently regular function -ij) on O, we will often employ the notations 
_ pij^] ( v-0)'^ = P'^[V'(/'] (V-i/' being R'^-valued, the projection is taken 

component per component). Further, for i<|L S £ and if*|L* G £*, we introduce 

(31) i/'a-ii, = /" V', V'a-il* = ^- / ^■ 

For L € cJQJt, there exists it|L C 951 that coincides with l; in this case we assign 
V'i = V'ail- If V'laa = 0, we have V'l, = for all l E dWl. For l* e we assign 

= J^^ ^. If ^ has a compact support in H, and size(T) is small enough, we 

have V'l-'^= for aU l* £ dTl* . 

Combining the above notation, we write -0^ = (]P[0], {''pK)KGdm, k*) k* £d<m*^ 
for the projection of a sufficiently regular function ip on the space R'^, and denote 
the corresponding projection operator by P'^. 
3.7. Dependency on t and further notation. 

• Let T be a DDFV mesh as described above. Let Ai > be the time discretiza- 
tion step. Set h — max{size(X), At}. By convention, we will use h as the parameter 
for a sequence of finite volume schemes; our interest lies in studying convergence of 
corresponding discrete solutions as ft, | 0. 

Denote by N the integer part of T/ At. In the sequel, in our notation we omit 
the dependency of iV, T and At on h. 

• For a functional space X on fi, we denote by S^* the projection operator 



-1 pnAt 

(32) : § € Li(0, T; X) ^ (S")„=i,...,Ar G {X)^ , §" = — / §{t) dt. 

J{n-l)At 
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• A discrete function on Q is a set it'^'^* = (u'^'")n=i,...,7v, where for each n, 
u'^'" is a discrete function on fl. The sej^of all such functions is denoted R^^'^. 

A discrete function on Q is a. set m'^''^* = {u'^'^)n=o,...,N7 where for each n, u^-" 
is a discrete function on fi. The set of all such functions is denoted by R,(^+i)x^. 
We also_use discrete functions u"^^^* e R^^^ and u^^^* G p^(w+i)xT_ ^jach of 
u"^'^*, u^'^^ ^u^'^* is therefore a restriction of u^'^*. The entries of u*^'" are denoted 
by (respectively, u^. ) for /f e 9JI U (respectively, for k* G QJl* U 9971*). 

A discrete field on Q is a set JT'^'^* = ^ where for each n, is 

a discrete field on Vl. The set of all such fields is denoted by (R'*)^'^^. 

• Any discrete function can be composed with a mapping A : R ^ R'" , 771 G N; 
for instance, A{u'^'^*-) stands for w'^"^* with values w"^ = A{u^) for k G 371 and 
w^, =^ A(u^. ) for K* G 971*, for n = 0, . . . , A^. Similarly, any discrete field can be 
composed with a mapping : R'' — > R™; one has v(-^^) = ('^(■^"))dgi>- 

• We say that a discrete function is nonnegative ( respectively, nonpositive), if 
all its entries are nonnegative (respectively, nonpositive); e.g., for v"^ G R'^ the 
notation v"^ > means that v,^ > for all k G SPT and v^* > for all k' G 97t*. 



3.8. The finite volume scheme. With the notation introduced above, the finite 
volume discretization of problem (1) takes the following compact form: 

find a discrete function m'^'^* on Q satisfying for n = 1, . . . , N the equations 



u 



x,n _ „,T,(n-l) 



(33) 



At 



+ (div cf)^ [w^'"] - div ^ [a( V^u;^^")] 



w 



A{u 



together with the boundary and initial conditions 



(34) 



for all n — 1, 




for aU K G dm 
for all K' G 9971*; 



(35) 
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Let US state (33) in a more explicit form: 



for all n = 1, . . . , iV, 

,("-!) 



—U 

m 



At 

d-1 



rrif 



size(T) 

K* San* 



n (1-1) 



m^-nK*(A(uK) -^(wa-)) = "1k§"-, for all k G aOJl, 



At 



. m^nK'iAiuj,*) - A{u^)) = m,,.§^.., for all k' e dWt* . 

SIZGI -i^ ) 



Here §" , S" are given by (30), (32); gK,L,gK*,L* are some numerical convection fuxes 
satisfying (28); i/k, i^k* for s given have the same meaning as in (24), (25); finally, 
for s given such that s C d G 2), V5A(it^'") is the vector of R'^ constructed from 
the values = ^(w^), w^- = MK') by formulas (21) (for d = 2) or (22) (for 
d = 3), i.e., in the way indicated in Remark 3.3. 

4. Elements of discrete calculus for DDFV schemes 

In this section, we list convenient formulations of various summation-by-parts 
formulas and chain rules needed for the analysis of the discrete problem (33). 
4.1. Discrete duality formulas for the diffusion terms. 

• Recall that R'^ is the space of all discrete functions on il. For m G N and 
w'^,v^ G (R^)", set 



1 d-1 



(36) w'^, v'^ ^ ;j X! ■ + — ~ ^ ' ■ ^'^^ 



- - d ^ a f-^ ^ 

(here • denotes the scalar product in R™); it is clear that |-, • is a scalar product 

on (R'^)™. We will use it for m = 1 or to = d. 

• Recah that (R'')^ is the space of all discrete fields on n. For ^'^,5'^ e (R'')^ , 

set 

(37) (.F^g^}}= ^m^J-^-e^; 

it is clear that "^'i is a scalar product on (R'')^. 

A key property of DDFV schemes (see [41, 11]) is the following discrete analogue 
of the duality between the — div [•] and the V[-] operators; it is sometimes called 
the discrete duality property for finite volumes. 
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Proposition 4.1. Let G RJ and T'^ e (R'')^. Then 

Proof. The proof is straightforward, using the summation-by-parts procedure. Let 
us give it for the case d — 2. Note that for d = s = s^i^., nis = ■^Tn^iLdKL = 
im^.i^.d,,.^.. By (36), by (24), (25), and finaUy by (21), (37), we get 

[-div^[.F^], ^ 



2 ^s<^e,,s=s1^ 



K'*G<Bt* sev*(K*) 



TOS >S • ( -j I^A'.L H !^A'*t* j 



'SG6 — 'nes LL 

Furthermore, we have the foUowing "entropy dissipation" inequaUties: 



□ 



Proposition 4.2. Let e RJ and -0 e I?(f^), > 0- iei 6i : R ^ R 6e a 

nondecreasing function. Assume that 

(38) either 9{Q) = 0, or t/' G P(r2) and size{%) is small enough. 

Denote tf^ = ¥^[i]]. Then 

div^[fc(V^A(ii^)) V^^Cu"^)] , 6l(w^)V;' 

< -|fc(VM(u^)) VMe(u^), V^^/r 



(39) 



Remark 4.1. Note that the conformity of the meshes (see Remark 3.2) is essential 
for this result, as well as the particular form of a and (for d — i) condition (19). 

Proof. Let us treat the left-hand side of (39) term by term. It is the sum of generic 
terms of the form T^^stT*, here 



T, 



K,s ^ 3 m^— m^|x.A:( Vs A(w'^)) Vs A{u'^) ■ Vk,lO{uk)^k 
a nriK 



1 



with s = s'^''" G v{k). The notation T*, ^ stands for analogous terms involving 
K* and s G v*{k'). Notice that thanks to assumption (38), 9{u'^)iIj'^ G RJ, so 
that we can also add the terms T^.s, T*, ^ corresponding to k G dWl, k' G 
respectively. The summation of T^.s, T*, g therefore runs on all subdiamonds s ~ 

s%'*,L* S ®i with the associated if, l G 371, if*, l* G 271*. 
The convexity argument yields 

(40) {A{z) - A{z))e{z) < Ae{z) - Ag{z) for aU z, z G R. 
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By (20), using the positivity of -0^ and applying inequality (40), we get 



(41) 



The terms Ti^-^s are treated in the same way. Now by the same computation as 
in the proof of Proposition 4.1, one shows that the right-hand sides of (41) and of 
the corresponding inequality for Tk*,s sum up to yield the right-hand side of (39). 
This concludes the proof. □ 

4.2. Summation formulas for the penalization terms. For the penalization 
operator , we have the following summation formulas. 



Lemma 4.1. Let w'^ G and ip"^ £ RJ. 

d-1 



The 



(42) 



E 



size(%) 



Further, Zef A, : R ^ R be nondecreasing. Assume u'^ E R^ is such that A{u'^) 
belongs to Rf . Let e V(Ti), tp > 0; denote = ¥^[il:]. Assume (38). Then 



(43) V[A{u^)],e{u^)4,^ 



d-l y-^ {A{uk)- A{uK')){i>K-i'K') 



size(T) 



In both formulas (42), (43), the values ip^ji^K* for K G StJT, if* G VCfl* are those of 
the corresponding discrete function ip"^ . 



The proof is straightforward from the definitions of 
summation-by-parts procedure. 

4.3. Discrete duality formulas for the evolution terms 



and 'P'^, using the 



Lemma 4.2. Let : R ^ R 6e a nondecreasing function, and ?7 = / d[s) ds be 

its primitive. Let tp G V(Q), tp > 0. Denote tp'^''^* = P'^ o S^^[ip]. Then for all 

^■xM g j^(7V+i)xx has 

N 



Y^At 



n=l 



,T,(n-l) 



^x,(n+l) _ ^,T,n 



At 



N-1 



> 



At 



Proof. The formula follows by the Abel transformation combined with the convexity 
inequality: {z — z)9{z) > r]{z) — rj{z) for all z, z G R. □ 

4.4. Discrete duality formulas for the convection terms. For the convection 
terms, we have a more involved "entropy dissipation" duality formula. For later 
use, we state it in the double framework, although each of the meshes 5(71,971* is 
treated separately in the proof. 
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Proposition 4.3. Let e RJ and ijj G T>{Vl). Let 6* : R ^ R 6e a nondecreasing 
function. Assume (38). Consider the associated entropy-flux pair 



rj^ J e{s) ds, q = e\- J f{s) de{s). 
Denote V'^ = P'^[V'] and ( V?/-)^ = P'^[V-0]. One has 
(divcf)"[u^], 0(u^)V'" 



(44) 



whe 



q(^.-), (v^)- +/,K-,^] + i?eK,v;] + /,*K',v] + i?;K*,V'], 



A'|i,e£ 



K*|l,*G£* 



with 



(46) 













J UK 


^9 






J u 



[gK.Lis, s) - gK,LiuK, u^)) dO{s), 

igi<',L'is, s) - gK'L'{uK',UL>)) dd{s). 



Further, one has Lg^'" > for all k\l G £, and the remainder term Rg satisfies 
(47) < (max|0(u,,)l) ^ m^^^ {R^^^- + Rf"-) {\^k-^kil\ + \^l~^kil\) , 



(48) i?^^ = |.gK-t(uK,UK) - gK-i,(w/f,Ui,)|, Rf'' = |.g/c,L(ui,, Wl) - .gK,i,(uK, Mi.)|- 

Similarly, one has Lg > for all A'*|l* G 6* , and the remainder term Rg 
satisfies the analogue of (47), (48) with K,L,Viyi,£ replaced by a:*, l*, f *. 

Note that our notation is consistent: we have /^'^ = /g'^, Rk'" = Rl'^ for all 
neighbours k,l (for dual neighbours k*,l*, similar identities hold). 

Proof. We exploit the ideas of [48] and [29]. 

Thanks to (38) and because is zero on boundary volumes, we have 

(49) 77(u^)V'l = for all L e aSUl; //(u^OV-l* = for all l' e 9071*. 

Separating the contributions of 9Jl and , we write the left-hand side of (44) 
as i/+ ^L*, where 

d a ' 



(50) 



Applying (28) (c) and (29), using (49) in the summation-by-parts procedure, we get 



K\L^£ 



/ = ^ ™A1L ( 6'(Ui,)(5A-,i,(Ui,,1ti) - ffir,i(uR',Mi)) ■01, 

-0{uK){gK,L{uK,UK) - gKL{uK,Ui^)) ijjj. 
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Hence, choosing i/'^ii as defined in (31), we have 

K\L££ ^ 

-^("K)(5K,i(wK,M/f) - g/f,I.(UA-,Wl,))^ 1pK\L 
-OiUK){gK,LiuK,Ui,) ~ gK,L{UK,U^)) {ip K - ''P K\l) )■ 



Now recall that g = gK,L satisfies (28) (b). Thus the following intcgration-by-parts 
formula holds true: 

(q(6)-q(a)).;.,,, = (0(6)f(6)-0(a)f(a)- f%{s) de[s)) ■ u^,, 

b 

= e{b){gib,b)-g{a,b)) - e{a){g{a,a)~g{a,b)) - / {g{s,s)~g{a,b))de{s). 

J a 

We deduce I = J + Ig + Re, where 

J ^ ^ TOif|t(q(wA-) - q(wt)) ■ Vk^i^iPk^l = ^ q(^tif) • ( X! "^ifii^J^fi^- 

51 ■ / 'J^^'<= J2 div(q(w«)?A) = J2 ■ '^''^ 

= ^ m^^- q(MK-) • ( V?/>)k, 

and 



lo ^ ^ m.K^L{ / (ffK,t(s,s) - gK,L{uK,u^))de{s)) i/jkil, 
< ^ m^-ii MgK,i,('«A-,UK) -5K,i,(uK,Mi)| + |5k,i,(ui,,Ui,) -gK,i,(MK,Wi,)| ] 



X + X (max |6I(uk)|). 

In the same way, /* = J* + Ig + Rg with analogous estimates. We have the equality 
i J + J* = q(u^) , ( V^Y . With the notation of (45)-(48) the resuh of the 
proposition follows. □ 



5. Properties of discrete operators and functional spaces 

In this section we state important embedding and compactness properties of 
spaces of discrete functions, as well as the asymptotic (as h ^ Q) properties of 
various discrete operators. 
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5.1. Discrete functions and fields as elements of Lebesgue spaces. 

For any E C Q, denote by ll^; its charaeteristic funetion. 
For 77, = 1, . . . , iV, set 

Qk = [{ri ^ ^)^t,nAt[xK, for G 

Qk' = [(" - l)At,nAt[xK% for K- e Tl*; 
Q'o = [{n - '^)^t,nAt[xD, for L> G X>. 

For a discrete function w'^'^* on Q, denote by w™-^* (respectively, by v™*'^*) the 
piecewise constant function 

N 

n=l KSajt 

("respectively, =E ^ llg. . (t, 

Whenever it is convenient, we identify the discrete function v"^'^^ G R^^'^ with 
the function on Q given by 

a a 

In a similar way, we identify a discrete field JT*^'^* g R^^^ on Q with the function 

N 

n=l DgS) 

Analogous conventions apply to time- independent discrete functions and discrete 
fields, in which case we suppress the superscript At in the notation. 

5.2. Consistency properties of discrete operators. In the proposition below 
we show the consistency properties of the projection and discrete gradient operators 
in Lebesgue spaces. Also note the property (w), which, combined with formula (42), 
expresses the fact that the penalization operator introduced in Section 3.4 vanishes 
(in an appropriate sense) as size(X) 0. 

Proposition 5.1. Let T be a double mesh of VI, At > 0, h ^ m.ax{size{'Z) , At} , 
and q E [1, +oo]. Then 

(i) there exists a constant C that only depends on VL, q and reg(T) such that 

Vw G L^iQ), ||(P" o§^*«;)™-^*||i, + ||(P^o§'^*«;)™*^'^*||i, < C||«;|U„ 

and 

Vu;Gi9(0,T;W(}'«(f)))), || V^P^ o <C\\\7w\\l^; 

(ii) for all w G L'^iQ), q < +oo, both {V^ o §^*'u;)'^'^'^* and (P'^ o §^*i(;)™*'^* 
converge to w in L'^{Q) as h — > 0; 

(ill) for all w G Li{0,T;W^'''{n)), q < +oo, the discrete fields V^P"^ o S^*w 
converge to Vw in {L'^{Q)Y as h —> 0; 
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(iv) let i) e V{n), and V-^'" P^(S[V'])", ?2 1, . . . , TV. There exists a constant C 
that only depends on Q and rcg(X) such that 

V At V m,,nK' ^'^"~t^'^' < C II WIIl- X sizeiT). 

Proof. The proof of (i)-(iii) is a straightforward gcncrahzation of [11, Lemma 3.3, 
Proposition 3.4 and Corollary 3.5]. We need to take into account the fact that 
IIS'^^'^llLnx) ^ IklU'CX) and (for q ^ +oo) ||§^*u; - w||l,(x) ^ as At 
for aU w G Li{0,T;X), where X stands for L9(S1) or for W^-'^{n). Remark 3.4 
is important for (iii) (thus, the Delaunay property of 3Jl is used). Further, in a 
standard way similar to [11, Lemma 3.3] one proves that for all k € k* € Wl* 
such that kCik' ^ 0, one has Itp^ - ^k* | < C(reg(T))|| Vt/'Hl^ x size(T) for all 
n — 1, . . . ,N. Hence the claim (iv) follows. □ 

5.3. Discrete embedding and compactness results. Next we state a version 
of the Poincare inequality and an cmbedding-kind translation estimate on double 
discrete functions. 

Proposition 5.2. Assume % is a double mesh on fi, At > 0. Let q G [1,+cx)). 
There exists a constant C > that only depends on diam (fi) and q such that 

(i) for all w^'^' £ R^""^ one has ||w^''^*||l. + ||u.™*'^*||l. < C \\ V^w^^^*||,; 

(ii) for all w"^ G Rq , for all A G R'' one has 

||^-(. + ^)_^™(.)||^, + ||^-*(. + A)-«;™*(-)|U, <C||V-z«-|U, X lAl^Z-^. 
Proof. The proof follows the lines of [12, Lemma 1] and [11, Lemma 3.6]. Note 
that if d = 3, the fact that all interfaces if]L are triangles plays an important role 
in the proof. □ 

Here is the asymptotic compactness result for "discrete i^(0, T; W^Q^'^(r2)" spaces. 

Proposition 5.3. Let p G (l,+oo). Assume we are given a family {w'^'^^^h of 
discrete functions in R^^'^ corresponding to a family of double meshes T such 
that reg(T) is uniformly bounded (recall that we parametrize the meshes by h = 
max{size(X), At} ). 

(i) Assume that there exists a constant C > such that 

Then there exists a (not labelled) sequence of meshes such that as h 

w'^''^* = — w™'^^ H w™* '^^ converge weakly in L^(Q) to some limit w; 

d d 

furthermore, w G £^(0, T; lyj '^(rj)) and 

the discrete fields V^w'^'^* converge weakly in {L^{Q)Y to S/w as /i — > 0. 

(ii) If, in addition, 



N 

E 

ri=l 



where are the penalization operators introduced in Section 3.4, then 

both w™'^* and w™' converge to w weakly in L^{Q) as h 0. 
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Remark 5.1. Note that upon providing uniform estimates on time translates of 



in LP((5), strong convergence to w in LP{Q) holds true (see Section 7). 



Proof, (i) The proof is very similar to the one of [11, Lemma 3.8]. 

First, by Proposition 5.2(i), both families {w'^' ''^*}h of components of 

are bounded in LP(Q). Therefore we can choose a common sequence such that 



w 



T.At _ 1 an, At 
d ^ 



-1 „,,aji*,At 



both components converge weakly in LP{Q). Also w 
converge weakly to some limit that we denote w. We can also assume that the 
corresponding sequence { V^i(j'^'^*};i converges weakly in {LP[Q)Y to some limit 
X. Let us show that w S Lp{0,T; W^'^in)) and x = Vt«. 

Take any field T G (/.^'(O, T; W^^-p' {^)Y . Denote by T"^^^^ the discrete field on 
Q with entries 



^ r-i 



1 



At X irLr 



riAt 



(n-l)At J D 



T. 



Denote by (divJF)^''^* the discrete function P'^ o S [div J^] on Q, which has the 
entries 



(div.F)^, 



1 



Atrrij. 



nAt 



K J {n-l)At J K 



1 



riAt 



T ■ V,, 



lAt 



AtrnK" J(„-l)At Jk 



By Proposition 4.1, by definitions of -Jj ■, 
in Section 5.1, we have 



1 



lAt 



AtmK' J(„-i)At ^^'^■^*(^'* 



and using the notation introduced 



TV 



TV 



ri=l 



= ^At|^^'" ' V^u;^'"| + ^At|div^[.F^'"] , 

n— 1 n— 1 

TV TV 

= E Ai{{^"'" , V-u;^^"}} + ^ A4(div^)-'" , , 

ra=l 

TV 

^At|div^[^^'"] - (div.F)^ 



(div^) + 



n=l 

■I At . Y7T,„x,At 



TV 



n=l 



^Ai div^[r^'"] - (div.F)-^'^ 



As in Proposition 5.1, one shows that — .F||^p' tends to zero as /i ^ 0. 

Therefore we deduce 



TV 



(51) 0= / J^-x+ / (divjr)u; + lim VAt div^[Jf^^^"]-(divJ^) 
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By definition of 

N 



n=l 



we liave 

^At|div -^[^-^'"] - (div^)^ 
1 

-"S^ AtirriKw'Z / (/ T - in„„T'^i) ■ Vk 

V " Aim, y(„_i)^,^sGvMV^^ 

VAi m,-.<. — / V (/ T^m„,T'^)-v^A. 



d-1 



Denote by i? + i?* the riglit-liand side above. Summing by parts, we get 

N ^nAt 



n=l k\l£s 



(n— 1) At J K\L 



T - 



1 



wliere d stands for the diamond containing the interface k|l. By the Holder 
inequality, we deduce that \R\ is controlled by 

N 



n=l fc|te£ 



At 

(n— 1) At J K\L 



1 



N 



"i, IP 



V \ — 



T„,,x.At 



At y ^ mK\LdKi. I — 

n— 1 K\L^e 

Using standard estimates similar to [11, Lemma 3.2] and the definition of V 
we conclude that 

\R\ < C(reg(X)) x size(T) x UHlp' (w^-p' }\\ ^^w^'^'\\l. 

< C(reg(T)) xhx miL.'iw^..') x C ^ 

as /i ^ 0. In the same way, wc find \R*\ ^ as /i ^ 0. 

Thus for aU T G {LP {0,T;W^-p' (n))'^ , the last term in (51) is zero, so that 
w e LP{0, T; W^'^i^)) and x = Vw. 

TV 



(ii) If also ^ At hP^[u;^^"],u'^'" 



< C, then by Lemma 4.1 with ip"^ ~ w"^ we get 



d-1 



TV 



^ At "^-n^* « - <. )' < t^^- 



This means that — w™*''^*||i2 — > as ft. — > 0, which permits to identify the 

weak limits of both ly™-'^* and ly™'^'^* with w. □ 

6. Properties of discrete solutions 
6.1. A priori estimates. 

Proposition 6.1. Assume we are given a family of double meshes T of D, and 
associated time steps At such that h = max{si2:e(T), At} 0. Assume that reg(X) 
is uniformly bounded. 

Let u^-^'^ be a solution to (33), (34), (35) (recall that w^-^* = A{u^'^'^)). Then 
the following a priori estimates hold uniformly in h: 
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(i) max{||u-^^*|U», < M := ||uoliLo= + / ||§(t, OIIlo^ dt; 

(a) there exists C > such that 



N 



V^w'^'^*||lp < C and ^Ai hP^[w-^^"],w^'" 



(in) there exists C > such that (with the notation of Proposition 4.3J 

AT 

^ At (/mK-'", 1] + /fdK"*'", 1]) < c- 

n=l 

(iv) there exists a modulus of continuity ij-'a(') such that for all A > 0, 
\w^''^\t + A, x) - w""^'^\t, x)\ + \w'^'^'^\t + A,x)- w'^''^\t, x)\ < uja{a), 



where w'^'^*' ,w'^* are extended by zero on {NAt,+oo) x il. 

Proof (i) Denote S* = (S'^* [§])'' and §^^* ^ P'^[(§'^* [§])']. For n = 0, . . . , iV, set 

c" = II^oIIl- +Er=i^i|l§iU~; note that c" < ||woI|l~ + \\Ht, -^l^ dt = M 
for alln = l,...,iV. 

Let us prove by induction that < c", < c". This claim is 

clear for n = 0. Assume it holds true for n = fc — 1. Take the scalar product 

of equations (33) corresponding to n = A: with the discrete function Oiu^'^ 
sign "''(u'^'''" — cJ^). We get 

- x,fc _ x,(fc-i) Ti rr — 

- - S^-^- , 6{u^^^) + (div,f)^[7.^^'^] , 6(u^^^) 



(52) 



At 



div^[a(V^A(w-^''^'))] , 6i(w^''^-) + 



,X.fe^ 



0. 



Let us apply to the last three terms above Proposition 4.3, Proposition 4.2 and 
Lemma 4.1 respectively, with ■0 = 1. Note that 0(0) = 0, so that (38) holds. We 
conclude that each of the three last terms in (52) is nonnegative. Hence 



> 



At 



> 
> 



At 
At 



S^Wl^-S^'") , sign + (u^'^--c'^-) 



sign + (w'^''' - c*-') 



)+ , V 



where 1^ = P'^[l]. By the induction hypothesis we deduce that {u'^''' — c'^)+ < 0, 
which proves our claim for n — k. 
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(ii) For n = 1, . . . ,N, take the scalar product 



of equations (33) with the 



discrete function w'^'" = A{u'^-"). Multiply by At and sum up in n. We get 



N 



At 



N 



(53) -^At |div^[a(V^u;^^")] 

n=l 
N 



n=l 
N 



, ^(u-'") +^At (divj)-^'"] , A{u-n 



+ 



n=l 



Note that with 9{-) ~ A{-) and ip = ^, (38) holds. Applying Lemma 4.2, Proposi- 
tion 4.3, Proposition 4.1 and Lemma 4.1, respectively, to the terms on the left-hand 
side of (53), we find 



(54) 



N 



r](u 



N 



Tl=l 



N 



N 



n—1 n—1 

,T,0\ 1 T 



W 



where Ba{z) = A{s) ds and Ia,Ia ^^'^ defined in Proposition 4.3. The first two 
terms in (54) are nonnegative; the next one is lower bounded by a constant times 



(II V 



T,,,X.At| 



ip) due to the coercivity assumption on a. By Holder's inequality, 



Proposition 5.1(i) and Proposition 5.2(i), the first term in the right-hand side of 
(54) is majorated by C(rcg(T)) x ||/||lp' x || S/'^w'^'^^Wlp- Finally, the last term in 
(54) is upper bounded by a constant times mo /_" ^(s) ds. Hence, (ii) follows. 

(iii) We proceed as in (ii), multiplying equations (33) by w"^'" instead of A(ii^^"). 
As in (54) above, taking 9 = Id, = 1, applying Proposition 4.2 instead of Propo- 
sition 4.1, neglecting the nonnegative terms on the left-hand side, we get 

N 



n=l 

N 

< ^ At |P^(§'^*[§])" , U^ "^ ^ 



n=l 



Using the L°° estimate (i) of the present proposition together with Proposition 5.1(i), 
we finally get (iii) with the constant 

C = C(reg(T)) xMx \\§U^ + ^m^ x (hoIlL-)'- 
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(iv) We adapt to the discrete framework the calculation that led to estimate (8) in 
the proof of Theorem 2.1(ii). Denote by J(a), J*(a), respectively, the integrals 

\u'"''^\t + A,x)~u""'^*{t,x)\ |A(u™^^*)(i + A,x)-A(M™'^*)(t,x)|, 

|u™*'^*(t + A,x)-M™*''^*(<,a;)| |A(M™*'^*)(i + A,x)-A(w™*^^*)(t,x)|. 
Let us first take fc G {1, . . . , N} and estimate the quantity 



N 



Jo(fc) := 



=k+l 



To do this, for n ~ (fc + 1), . . . , iV we take the sum in i from (n — fc + 1) to n of 
equations (33) and make the scalar product ■ , • with the discrete functions w'^'", 

where u'^'" := A{u'^'") - G RJ for n ~ {k + 1), . . . , N . Summing in n 

and assigning w'^'" = for n = 1, . . . , A; and n = (N + 1), . . . , {N + fc — 1), we get 

N 



Jo(fc) 

At 



E 

n=k+l 

N n 

e 



(55) 



n—k-\-l i—n — k-\-l 

N min{fc.N-?.+ l} 

E E 

1=2 j=max{l,fc-j+2} 
fc N 



,T/j _ ,,x,(i-l) 



At 

At 



At 



V 



= ^^At|- (div,f)^[u^-*] +div^[a(V^u;^'0] 
j=i i=i 

We claim that the right-hand side of (55) is bounded by a constant independent of 
h. Indeed, for each j = 1, . . . , fc, define Zj"^'* = , i = 1, . . . , iV. First, from 

the property (ii) of the present proposition and from formula (42) we deduce 



N 



(56) iiv^zj'^*iu. <c, E^Hr't^ri 



j=l 



< C, for aU j = 1, . . . , fc. 



In the sequel, we will omit the dependency of the entries of Zj'^^ on j. 

By definition of (div cf)'^ [■], taking into account that zj'" e Rq and using 
summation-by-parts, we deduce that for all j = 1, . . . , k, 

N 



^3 



^At|(div,f)-[^.-'^], z; 

i=l 

^ /I 

E'^M^ E '«ail3k-i,«,<) (4-4) 

XI «.,<.) (4.-2^)j 
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Since by (i), u™.^*^ u™*^'^* are bounded by M, using property (28)(d) we bound 
all values of gK\LTgK*\L' above by Cu!m{M). It follows by Remark 3.3 that 



where S = S^^l Hence 



N 



■h,j < C(d- l)tJM(A/)^At^ msl Vsz'^'''| < const || V^zJ'^^'Uli. 



i=l sec 



Using (56), we can uniformly bound Jij. Further, by Proposition 4.1 and the 
Holder inequality. 



J2,j ■- 



N 



^At|div^[a(V^w^'*)] , z 

1=1 



<||a(W^*)IU,.l|VX'^*l 



LP- 



Using the growth assumption on a together with (56) and (ii) of the present lemma, 
we can uniformly bound J2,j- Next, by (42) and the Cauchy-Schwarz inequality, 



J: 



■S.J 



N 



i=l 



N 



d— 1 \ ^ X ^ 



< 



d 

d-l 



i—1 K^DJl.K* GOJl* 

N 



^size(T) ^size(T) 

1/2 



1=1 



N 



size(X) 



1/2 



X 



E-^ E 



N 



k'Got.k'* Got' 
\ 1/2 . ^ 



size(X) 



1/2 



T,4 

z ■ 



Using again (56) and (ii) of the present proposition, we can uniformly bound Jz.j. 
Finally, like in (54), we have 



N 



1=1 



EAt p"(§^*[s]r , 



<C(reg(T))x|lS||^,. x||VX'^*|U. 



which is also uniformly bounded, thanks to (56). Gathering the estimates above, 
we conclude 

Jo(fc) < At E ( -^iJ + "^2 J + Jsj + ^4 J ) < C kAt. 



36 



B. ANDREIANOV, M. BENDAHMANE, AND K. H. KARLSEN 



Using the definition of 



and the L°° estimate on m™'^*, cf. (i), we get 



1 d-1 f^^^ 

(57) -J{kAt) + ——J*{kAt) < Jaik)+ / m^M max{±yl(±il/)} < C kAt. 

d " J(N-k)At 

Now let < A < r. We have a/ At = {k - I) + a for some k G {1, . . . , N} and 
a g [0, 1). Since u^'--^* is piecewise constant in t with step A^, we have 

J(a) = J((fc - l)At + aAt) < aJ{kAt) + (1 - a)J{{k - l)At) 

< aC kAt + (1 - a)C {k - l)At < C {{k - 1) + a)At = C A. 

From (58); together with the calculation used to pass from (8) to (9) (cf. the proof 
of Theorem 2.1), we deduce the required estimate 

\A{u'^'^'){t + A,x) ^ A{u^'^*){t,x)\ <t^A(A), A>0. 



(58) 



□ 



Q 

Similarly, time translates of A(m™*'^*) are controlled with J*{kAt) in (57). 
6.2. Existence of discrete solutions. 



Proposition 6.2. Let % be a double mesh offl and At > 0. There exists a solution 
yi,At of the finite volume scheme (33), (34), (35). 

Proof. First note that it is sufficient to prove existence of solutions wj''^* to (33), 
(34), (35) with A[-) replaced by a strictly increasing function Ap{-). Indeed, using 
the L°° estimate (i) of Proposition 6.1, which is independent of the choice of A(-), 
we get compactness of w^'^* in the ffiiite-dimensional space p{,(^+i)x'5^. Choosing 
a sequence of strictly increasing functions Ap that converges to A uniformly on all 
compact of R, we pass to the limit in the scheme (33), (34), (35) written for (a 
subsequence of) Ap{-) and Wp'^* and obtain existence for general A{-). 

Let us now assume that A{-) is invertible and rewrite the scheme in terms of 
with w'^'^* = A~^[w'^'^*). The existence of w'^'" is shown by induction on 
n = . . . ,N . For n = 0, solution is given by (35). Assume that exists. 
Choose • , • as the scalar product on R'^. We are looking for a solution tu^'" to 
= 0, where the operator L is given by 

A-i(z^) - A-i(u;^^("-i)) 



L : e R' 



At 



+ (divJ)-[^-i(z-)] 



-div^[a(V^2^)]+7'^[z-^] -P-^(§^*[§]) 



By Proposition 4.3 with 6* = Id and = 1, by Proposition 4.1 and by Lemma 4.1, 
there exists a constant C = C( 



— 1 1 



RX,||P-(S^*[§])"llR-,Ai such that 



By the coercivity assumption on o and by Proposition 5.2(i) we have 



(59) 



a( V^z^) , V^z^ \ > const || V^z'^H^^ > const 



\LP) 



Because the right-hand side of (59) is equivalent to (||z'^||r^t)'', we conclude that 
L[z'^] , z'^ > for ||z'^||jix sufficiently large. The existence of w'^'" follows by 
the standard Brouwer fixed point argument (see [G4, Lemme 4.3]). □ 
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We point out that the uniqueness and, more generaUy, continuous dependency 
of the discrete solutions on the data can be estabhshed as well (see [48, 49, 11] for 
results of that sort). However, in view of the convergence result of Theorem 7.1 
and the wcU-poscdness of the continuous problem, we view these questions to be of 
less importance. 



6.3. Discrete entropy inequalities. 

Proposition 6.3. Let % be a double mesh of and At > 0. Consider a solution 
u^'^* to the scheme (33), (34), (35); recall that w"^'^^ = A{u^-^^). 

Letip G V(Q), ip > 0; seti)^-^^ = W o'§>^^[il}]. LetO -.R^R be a nondecreasing 
function; assume that and 9 are chosen so that (38) holds; assume that At is small 
enough. Then 



N-l 



N 



T,(n+1) _ ^T,n 



At 



N 



(60) 



n=l 



^At|q(w^'"), (VV)'^'"J -^At|fc(V^w^^") V^Ae(w^^"), V^V^'" 

1 ri=l 
N 

(""'°) , ^"^i +EAtfp-(s^*[/])", 0(u-^")^-' 



N 



d — 1 \ ^ X ^ 

> — — 2^ At ^ m^r 



«-<,)(C-^^) 
size{%) 



where A0{-)^ A0{-) and ri{-) , c\{-) , Rg[- , R*g[- , ■] are introduced in Definition 2.1 and 
in Proposition 4.3, respectively. 

Moreover, with the specific choice 9 = 1 and tp G I?([0,T) x 17), there holds 



N-l 

n=l 
N 



^T,(« + l)_^l 



N 



At 



N 



(61) 



^At|fc(V^«;^'") V^w;^'", V^7/.^'"}> + ^At p-^(§^*[§])", 



N 



d-ls^. ^ «-<.)(V'5J-V'^) 



n=l 



K^Wi.K* eon* 



size{%) 



Finally, with the specific choices 9 = A and ip = C(i)j where e I?([0,T)) 
is a nonnegative, nonincreasing function with (^(t) = 1 for small t, we have with 
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B{z) = J^A{s)ds 



(62) 



N-l 



^T,(n+l) ^T.n 

At 



n=l 

N 



Proof. Inequality (60) follows by an application of Lemma 4.2, Proposition 4.3, 
Proposition 4.2 and Lemma 4.1. Note that in (60), we have neglected the positive 
terms Ie[u^-'", ■0"], Ig[u^*'", ip"']- In (61) the corresponding terms are zero because 
9 = 1, and we use the equality of Proposition 4.1 instead of the inequality of 
Proposition 4.2. Also notice that the term with ip'^-'^ in Lemma 4.2 disappears 
because At is small and i/j vanishes in a neighborhood of t = T. Finally, in (62) we 
have treated A{v^'^*) ^'^^ as a mere test function by applying Proposition 4.1 on 
the right-hand side, but we have used Lemma 4.2, Proposition 4.3 and the choice 
of the constant in x function ip"^'^*- to deal with the remaining terms. □ 

6.4. Control of the remainder terms in Proposition 6.3. For all £ T^{Q)i 
the terms on the right-hand side of (60), (61) coming from the penalization operator 
vanish as ft, ^ 0. Indeed, using the estimates of Proposition 6.1(i),(ii), the Cauchy- 
Schwarz inequality. Proposition 5.1(iv), and the boundedncss of 9 on [— M, M], we 
obtain 



N 

IE-* _E_ 

/ N 



9{u^) 



size(T) 




size(T) 



1/2 



< C|| Vi/'IIl- X size(T). 

Let us show that the terms ^^[w™* , "0] in (60), (61) (which are defined 

in Proposition 4.3) vanish as ft ^ 0. This holds true thanks to their upper bounds 

in terms the quantities /id[u™,l], quantities which are controlled by 
means of Proposition 6.1(iii) (known as the "weak BV estimate", cf. [48]). 

Proposition 6.4. Let g^^L G C(R^) be a function with properties {28)(a),(d). For 
a, 6 G R, consider 



jK\L 



ia,b) 



{9k,l{s, s) - gK,L{a, b)) ds, 



RT{a,b) = \g^^r.ia,a) - gK^r.{a,b)l R'^^{a,b) = \g^^^{b,b) - g^^^{a,b)\. 
There exists a continuous strictly increasing convex function Hm '■ R-^ R^ 



that only depends on C and ujm{') in {28)(d) such that 11^/(0) 
the following bounds hold: 

r RT{a,b)<Ul^{I,f{a,b)), 

\ i??-(a,&)<n^/(/j-^^(a,6)). 



0, U\j{0) = and 



(63) 



for all a, be [~M,M] 
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The proof is based upon the foUowmg generalization of [4N. Lemma 4.5]. 

Lemma 6.1. Let g G C{[a,b]) be a nondecreasing function equipped with a modulus 
of continuity to. Then 

{g{s)~g{a))ds> / cj'Hr)dr. 



Proof. Set S = ll! ^{g{b)—g{a)). Since \g{b)—g(s)\ < uj{b—s) and g is nondecreasing 
we have 

5(6) -cj(fe-s), b-6<s<b 
g{a), a < s < b — 6. 

Hence setting z — b — s, integrating by parts, and setting r = lu{z), we deduce 



.9(s) > 



(gis) - 9{a)) ds > / {g{b) - g{a) ~ uj{b ~ s)) ds 

Jb-5 

= 6u!{S) ~ / uj{z)dz= / zdijj{z)= / u!^^{r)dr. 

Jo Jo Jo '-' 

Proof of Proposition 6.4. Consider the case a < b. By (28)(a), we have 

Iidi^i ^) = / igK.Lis, s) - gjci(a, 6)) ds> {gK,L{s, b) - gnA^-^ ^)) ds; 
applying Lemma 6.1 to (/(•) = gK,L{',b) and recalling (28)(d), we deduce 

fgK,Lib,b)-gK,Lia,,b) p r'^^'" {a ,b) 

hTia,b)> {CuJM)-\r)dr^ {C uo M)-\r) dr . 

Jo Jo 

Thus in order to estimate R^'^{a,b) as in (63), it is sufficient to take the function 

nR 

III/ : i? G / [Clom) ^(r) dr. Clearly, IIm is continuous, strictly increasing, 

convex, HmCO) = o" and n^/(0) = 0. 

The other estimate in (63) is obtained in the same way, and the case a > 6 is 
obtained by symmetry. □ 

Corollary 6.1. (i) Consider Iid[u™ ,1] defined as in (45), (46) with 9 = Id, and 
■0 = 1. For general nondecreasing 0{-) and general G T^(S^), consider i?6i[M™,V'] 
defined in (47), (48). Assume \\u™\\oo < M. Let IIm be the function given in 
Proposition 6.4. LetT\\j be the conjugate convex function of Tim . Then 

(64) 

< m\ca-M,M]) inf f ^^/mK\ 1] + -n^,(2a max ^-^tlptu 

where C depends on reg(T), d and Q. 

(ii) Assume we are given a sequence of meshes T with size{%) — > and time steps 
At 0. Letu"^'^^ be the corresponding discrete functions such that ||u™''^*||oo < M 
and 'Ylin=i^t 1] < C uniformly in T, At. Choose V' G T^{Q) cmd take 

^" = (§'^*[,/;])". Then At R*g[u^''\ V"] ^ as sizeil) 0. 

Analogous statements that involve ^^=1'^^ -^wl*^™*'!] '^'^'^ '4'k-', '4^k*\l* with 
K' em*, L' eM*{K') hold for Y.'^^^ At 
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Proof, (i) By (47) and Proposition 6.4, for all a > we have 



< 2||6i||c([-M,M]) (^^^/fii-d/cAii^ X nj,/(/j'; 



dK,K\L 



Note that d^L < sizc(T). Further, even in the case the diamonds are not necessarily 
convex, the definition of reg(X) permits to control the multiplicity of the covering 
of by the convex envelopes of k and i^, k e 9Jl, l e m{k). Thus one can upper 
bound Y^Keon, l&u(k) 'mK\LdK,K\L by C(reg(X), d) mn- Applying the inequality rs < 
nAf(r) + n^^(s) on the right-hand side above, we deduce (64). 

(ii) First notice that for all G T>{Q), there exists C > such that 

max ' ^ < C, for all h > 0. 

n = l,...,N, K£SJt,L£^r{K) dfc.KlL 

Applying (i) for each n and summing over n = 1, . . . , A^, we get 

N / . N 



(65) 



y At V"] < C inf f y At /id[u™, 1] + T-UlJCa) 

^-^ a>o \ a ^-^ a 

n=l \ n=l J 



Q>o \ a a ^ 



where C stands for a generic constant independent of h. 
We have (nj,/)'(0) = 0. Therefore 

(nl,)'(0) = hm inf (a < lim ( b ^) = 0. 

b^O a \ b J 6-.0 \ b J 

Hence for all C > 0, lima^o ^njj(CQ;) ~ 0. We deduce that the right-hand side of 
(65) tends to zero as size(T) — > 0. □ 



Remark 6.1. Notice that if f is locally Lipschitz continuous, both Hm and H^j 
are quadratic; thus we can bound by Const ft.'' for all (3 < 1/2. Using 

the Holder inequality instead of the Young inequality, one recovers the result of 
[48] with f3 = 1/2. Whenever f is locally Holder continuous of order 7 < 1, we find 
n*,f(s) = Const fii+T. It follows that \Rg[u^,^lj]\ < Const /j/^ with /3 = under 
the assumptions of Corollary 6.1(ii). 

6.5. Approximate continuous entropy inequalities. Relying on Proposition 
6.3, we now deduce the limiting (as h —^ 0) entropy inequalities and the limiting 
weak formulation; one should notice that they continue to hold if we replace (77^ , q^) 
by regular "boundary" entropy-entropy flux pairs {rj^^, qj^). 

Proposition 6.5. Consider a family of double meshes % of and associated time 
steps At > 0, parametrized by h = max{si2:e(T), At}, h ~> 0. Assume that reg(T) is 
uniformly bounded. Denote the corresponding discrete solutions of (33), (34), (35) 
by u^-^K Fix € V{[0,T) xU), ip > 0, and set V''^-^* = P"^ o §^*[-0]. Fix as 
one of the functions rj^ , c S R. Assume either {c,^) £ R* x P([0, T) x ft), or 
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(c,^) eRxV{[0,T) X n). Then 



liminf , , , 
h^o V Jq d 



+ / ^(q^(^™^^*) + id-l) qtiu"^'^^')) ■ VV' 
(66) - / ki\7^w^-^')V^ A^^±yiw^'^') ■ W 

•J Q 

^ IJd 0'^^ + - 1) ^c^(^™'^°)) ^-(0, •) 

+ ^ ^(('?c^)'("™''^*) + ('^ - 1) ('yc^)'(«™*^^*))§'A) > 0. 

Furthermore, if € P([0,T) x il), we have 
(67) 

+ y" Q(f(w™'^*) + (rf-l)f(w™*^^*)) -fc(V^u;^'^*)V^w^^'^*j • 

+ y i(u™-o+(d-i)w™*'°)7/.(o,-)+ y" sv^)=o. 

Proof. By the choice of (c, ?/>), (38) holds. Thus, by Proposition 6.3, (60) and (61) 
hold; it suffices to develop these formulas using the definitions of • , • ' l^' ' 'J'- 

The second term in (60) rewrites exactly as the corresponding term in (66). 
Regarding the other terms on the left-hand side, we also use the uniform bound 
on m'^'" in L°°, the uniform bound on A;( V^w;^'^*) V^w^'^' in {Lp'{Q)Y, and the 
convergences 

E E ^-"^^ V. - O^i., E E J^-" ^Q".. - d^i. in L^(Q), 

n— 1 ftT^an n— 1 /cGot* 

TV N 

^ ^(P-(S^*[§])").^11q.-S, ^ ^ (P-(S^*[S])'%.llQ.,--SinXi(Q), 

n=l Kesnt ri=lA'*G53i* 
TV AT 

EE^JJ^Q^^^'E E iig^. -^ini°°(g), 

yT^T,At ^ ^p(g) ^-(0, •), V'™''^(-) V'(0, •) in L^l^), 

as /i — > (here we have put Proposition 5.1 to use). Finally, the terms on the 
right-hand side of (60) vanish as — > 0, thanks to the initial remarks made in 
Subsection 6.4 and Corollary 6.1(ii). In the same way, (67) follows from (61). □ 



7. Convergence and statement of main result 
We are now in a position to state and prove the main result of this paper. 
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Theorem 7.1. Consider a family of double meshes Tof^l and associated time steps 
At > 0, parametrized by h ~ max{sz2e(X), At}, ft, — > 0. Assume that reg(T) is uni- 
formly bounded. Then the corresponding discrete solutions m'^'^* of (33), (34), (35) 
exist, are uniformly bounded, and converge to the unique entropy solution u of (1) 
in the following strong sense: 



u 



m.At , „, „,OT*.Ai 



u, u in L'^iQ) for any s < oo, 



V^w'^'^* -> Vw in LP{Q), where w = A{u). 

Proof. Wc follow step by step the proof of Theorem 2.1. 

(i) Discrete solutions m*^'^* exist by Proposition 6.2. Besides, they verify the 
asymptotic entropy inequalities (66) (where we can replace rj^ by rj^^) and the 
asymptotic weak formulation (67), both of Proposition 6.5. 

(ii) Proposition 6.1 yields uniform estimates on both u™'^^ and u™''^^ in 
L°°(Q); on the time translates of both w™'^^ and w™' in L^{Q); on the pe- 
nalization term Y.n=i w'^'" ; and on V^w^'"^* in LP{Q). The latter 
estimate implies further uniform estimates: namely, an estimate of the space trans- 
lates of both u'™''^* and lu™*'^* in L^{Q), by Proposition 5.2 (ii); an estimate of 
V'^yl^^i y{w'^'^*) in LP{Q), because A^^±^y{-) is Lipschitz and by construction of 

V'^[-]; and finally an estimate of a( V^iw'^'^*) in (Q), because of the growth 
assumption on o. 

(iii) Thanks to the estimates of (ii), there exists a (not labelled) sequence of 
T,At with /i — !■ such that 

• by the Frechet-Kolmogorov theorem, each of the sequences w™'^* and w™* 
converges strongly in L^{Q) and pointwise a.e. in Q; 

• by Proposition 5.3, the limits of w™'^* and ui™*'"^* coincide (we denote the 
limit of ly'^'^*,?!'™*'^* by w), and V^w'^''^* converges weakly in LP{Q) to Vw; 

• a( V^u''^''^*) converges weakly in LP {Q) to a limit field x; 

• the sequences u™'^*, u™*''^* converge to : Q x (0, 1) — > R, respectively, 
in the sense of nonlinear L°° weak-* convergence (5). Also by (5), the functions 
yjm,At _ ^^,^!Bt,At) converge to A{ii) in the L°° weak-* sense; since the functions 
yj<m,At ^Yao converge strongly, A{pL) is independent of a and coincides with w. In 
the same way, we deduce that A{^*) is independent of a and coincides with w. 
Also observe that u™*'" both converge to uq a.e. in and in i^(ri). 

(iv) As in the proof of Theorem 2.1, we can use the chain rule and the Green- 
Gauss formula to deduce 

'-,{m + {d-mi^*))-^A{u) 

d 

''ifi^i)■VAi^) + {d-m^i*)■VAi^i*))=^ [ A^iw)-n = 0, 

Jo Jon 



(68) ^ 



Q 
d 



where A^ is defined in (15). 

(v) Next we pass to the limit in (67). Indeed, by (iii), 



1 1 



(gg) , dtu + div - (f(AO + {d- da = divx + S 



in LP {0,T-W-^'P' {^)) + L\Q), 



u\t=o = uq- 
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where 



u{t,x) ~ j fi{t,x,a)da, /i = — (/i + (d — 1)//*) . 



(71) 



Let us identify x (the weak limit of o( V^w'^^'^*)) with a( Vw), and consequently 
obtain that the weak convergence is in fact strong in U'{Q). To this end, we will 
use (iv) and (62) to establish the inequality 

(70) / X- Vio > liminf ^At |a(V^«;^'"), V^w"^'"!- 

•'Q ~^ n=l 

Indeed, using (69), we can represent the left-hand side of (70) as 
X • Vw C = - / {dtu, w C) 

+ I' (f(M) + id^ l)f(M*)) -VwC + J^SwC, 

where C £ 2?([0,T)) is nonincrcasing with C(t) = 1 for t small. 

Note that since A is nondecreasing, since u(t, x) is a convex combination of the 
values ^i{t,x, ■) and fj,*{t,x, •), and because A{^) = A{fj,*) = w, we conclude that 

w = A{u). 

To control J^{dtu,w0, we argue along the lines of the proof of Theorem 2.1. 
The duality product {dtU, A{u)) is treated via the weak chain rule (cf. [4]). Hence, 
exploiting also the convexity of B{z) = A{s) ds, 

'^{dtu^Aiu) 



(72) 



B{u)C' - / Biuo) 



fi{t, X, a) da ) (— C') 



< - / C' / B{'jl{t,x,a))da 
Jq Jo 

Using (68) and (72), we deduce from (71) that 



Biuo) 
Biuo). 



X- VwC> / C BiJlit,x,a))da+ / S(uo) + 



(73) 

On the other hand, Proposition 6.3 permits to evaluate the right-hand side of 
(70) as follows: 

/ N _ _ _ \ 

lim inf XI V^w^'") , w^'" C^'" ' 



(74) 



<liminf(EAt|B(.--"), ^ 



Biu^'") , r- 



N 



n=l 



X At P^(§^*[S])", W^'"C' 
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By the previously established convergences (see also the proof of Proposition 6.5), 
the right-hand side of (74) is equal to the right-hand side of (73). Once we let ( 
tend to ll[o,T)5 this establishes (70). 

Starting from (70), we apply the Minty-Browder argument that we employed for 
the continuous problem in the proof of Theorem 2.1. 

Take v e Lp(0, T; Wg •^(f])) n L°°(Q), and set v^'^* = o §'^*[z;]. In view 
of (70), taking into account the strong convergence of V'^w'^'^* to Vv in L^^Q), 
cf. Proposition 5.1, and the monotonicity of a(-) we obtain 

N 



I X - V(u'-t>) > liminf ^ At |a(V^w-^^"), V^u;'^^" - V^w^-"| 

■'Q n=l 

TV 



As is a well-known property of Leray-Lions operators, the strong convergence of 
Y'^^T.At L^{Q) implies the strong convergence of a( V^w"^'^*) to a( Vw) 

in LP'{Q). Therefore (75) yields 



/ X ■ V(w -v)> / a( Vw) • V{w - v). 
Jo Jo 



IQ JQ 

Choosing v = w ± X^j with A j and ?/> e LP{0, T; W^'^'ifl)), we conclude 

X = a( Vu;). 

Moreover, as in the proof of Theorem 2.1 and [11, Theorem 5.1], relying on the 
strict monotonicity of a and utilizing an argument of [24, 21], we also deduce the 
strong convergence of V'^u>'^'^* to Vw in U'[Q). 

(vi) Now we can pass to the limit in the weak and entropy formulations listed 
in Proposition 6.5. The passage from (67) to (D'.2) is straightforward. In (66), 
we first work with regularized boundary entropies. Taking the limit, all the terms 
converge to the corresponding terms in (D'.3) in a straightforward way, except for 
the third one. Let us show that 

V'^A^^±^^,(w^''^*) converges weakly to VA(^±^j,(w) in L'P{Q). 

Indeed, both A^^± y(u>™''^*) and A^^± y{w™''^*) converge to ^(,,± yi'w) by the 

a.e. convergence of w™'^* ,w™' to w and the continuity of ^(^i y- Using the 

boundedncss in U'{Q) of V"^ A^^± y{w'^'^*)) and the compactness property of 
Proposition 5.3, we conclude that our claim holds. The subsequent arguments 
arc the same as in the proof of Theorem 2.1. 

(vii) We conclude that (^,/i*,w) is an entropy double-process solution of (1). 
In view of Theorem 2.2, this brings to an end the proof of Theorem 7.1; indeed, 
we obtain the convergence to u for each sequence of discrete solutions with — *■ 0. 
Also, the fact that ^ and /i* turn out to be independent of a means that the 
convergence of u™'"^*, u™'^* to u is strong in L^{Q) for all finite s. □ 



8. On the choice of FV scheme and various generalizations 
In this section we discuss other possible choices of finite volume schemes for (1). 
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• The use of DDFV schemes is motivated by their convenience when it comes to 
the discretization of nonhnear diffusion operators. Other possibihties exist; among 
them, let us mention the schemes studied in [56] (sec also [3, f 2, 8]), in [10], in [42] 
(see also [47]), and in [52] (see also [50, 51, 53]). All these schemes possess some 
variant of the "integration- by-parts" property of Proposition 4.1. 

The 2D schemes of [10] are restricted to Cartesian meshes, so they do not allow 
for domains much more general than rectangles. Notice that their generalization to 
3D appears to be straightforward. The techniques used in the present paper and 
in the references we cite, such as [48, 49], combined with those of [10], allow to 
design rather simple FV schemes on Cartesian meshes for problem (1) and to prove 
their convergence. In this case, the notion of entropy double-process solution is not 
needed, and the theoretical results in [40] can be adapted directly. 

This is also the case of the "complementary volumes" schemes as described in 
[56]. In 2D, ideas quite similar to that of [5(i] were used to construct the schemes 
of [3, 12, 8]. All these schemes work on meshes dual to conformal triangular 2D 
meshes, and the discrete gradient is reconstituted by afhnc per triangle interpola- 
tion. "Complementary volumes" schemes are simpler than our DDFV scheme from 
the practical point of view, since one discrctizcs the problem on the same mesh 
SPI* using, roughly speaking, half of the unknowns. The discrete duality properties 
for the 2D "complementary volumes" scheme are shown in the same way as for 
our DDFV schemes; the proof is based upon Lemma 9.6 (see Appendix B and also 
[7, 8]). Unfortunately, the straightforward generalization of these "complementary 
volumes" schemes to 3D fails to satisfy the discrete duality property, except for 
very constrained geometries of the meshes (see Remark 9.2). 

The key feature of the 2D schemes of [3, 56, 12, 10, 8] (see also [40]) lies in the fact 
that the fluxes across interfaces are reconstructed "manually" . The approaches of 
Droniou and Eymard [47, 42] and those of the HVF, SUCCES and SUSHI schemes of 
Eymard, Gallouet and Herbin [50, 51, 53, 52] arc different; they rely on introducing 
additional unknowns (either for the fluxes, or for the values on some of the edges) 
and on careful penalization of the finite differences. 

The schemes HVF, SUCCES and SUSHI (among many others) were designed for 
handling linear anisotropic, heterogeneous diffusion problems with possibly discon- 
tinuous coefficients; in this framework, their convergence is justified. These schemes 
avoid usage of double meshes and thus may have less unknowns; they work both in 
2D and 3D. We refer to Eymard, Gallouet and Herbin [52] for the description and 
comparison of these and related (e.g., mimetic finite difference) schemes. Finally, 
let us also mention the schemes of Aavatsmark et al. (see, e.g., [1, 2]), that are 
in a sense intermediate. The gradient reconstruction used in [1, 2] also involves 
additional edge unknowns, which are eliminated by solving, locally, an algebraic 
system of equations. 

The scheme of [42] designed for nonlinear Leray-Lions kind problems can be 
directly compared to the DDFV schemes of [11] and of the present paper. The 
scheme of [42] is very interesting because of the extreme generality of the geometries 
allowed for the mesh (and it works in any space dimension). For this same reason, 
theoretical justification of its convergence in the hyperbolic-parabolic framework 
(1) seems problematic. Indeed, the conformity (orthogonality) condition was used 
in an essential way in the derivation of the discrete entropy inequalities (see Remark 
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4.1). The same difficulty arises for tire double schemes of [ ] in the case of non- 
conformal meshes, cf. Remark 3.2. In passing, let us point out that the conformity 
(orthogonality) assumption on the meshes is the only condition that is known to 
ensure the discrete maximum principle for the DDFV schemes. 

In conclusion, the 2D and 3D conformal DDFV schemes studied herein, although 
constrained by the orthogonality condition, by the Delaunay condition, and by 
condition (19), combine some degree of flexibility (e.g., any polygonal/polyhedral 
domain can be partitioned into triangles/tetrahedra satisfying these restrictions) 
with the rigid structure properties underlying our convergence proof. But because 
of the conformity constraint, the advantage of simple local refinement procedures 
for 2D DDFV schemes, pointed out in [11], is lost. 

• Our assumption that 9JI consists of simplcxes is a practical one simplifying 
the presentation of the scheme. In 2D, it can be replaced by the more general 
assumption that any element of £CJt admits a circumscribed circle. In 3D, we can 
assume that each k G QJl admits a circumscribed ball, and each interface k|l is a 
triangle satisfying (19). 

Notice that Remark 9.3 (see also [7]) makes it possible to define a consistent 
discrete duality scheme even when the interfaces k|l are not necessarily triangles. 
Unfortunately, the discrete Poincare inequality may fail in this generality; this 
undermines the subsequent convergence analysis. Yet one interesting case is that 
of a Cartesian mesh QJl; the corresponding DDFV schemes are alternatives to the 
scheme of [10] discussed above. More generally, one can start with a mesh 9Jl made 
of rectangles (e.g., inside fi) and triangles (e.g., near the boundary dO) in 2D. 

• As pointed out in Remark 3.5, a different kind of reconstruction formula is 
needed for problems in 4D and higher dimensions. It would be interesting to con- 
ceive discrete gradients consistent with affine functions, following the principle for- 
mulated in Remark 3.3. One natural way is indicated in [3.s]. 

• The choice of penalization in our double scheme can be changed (sec Remark 
3.6). One could also penalize the differences {uk — Uk*) instead of the differences 
{wk — Wk*) = {■A{uk) — ^(uft-)); this would permit to avoid the use of double- 
process solutions. But this choice would introduce additional coupling between the 
sets of variables {uK)Kemt and {uK*)K*eoyi* in the "hyperbolic" regions. Indeed, 
if, e.g., A{u) = 0, there is no coupling at all between the variables sitting on SOT 
and those sitting on SUT*. Therefore our choice seems more convenient in terms of 
practical implementation. 

• Convection-diffusion problems with anisotropic linear and nonlinear diffusion 
were considered in [35, 34] and in [17, IH]. General DDFV schemes do not seem 
easy to adapt to the nonlinear anisotropic framework, because of the presence of 
"privileged" directions of diffusion. In this case, the schemes of [10] on Cartesian 
meshes constitute a natural choice, and the geometry of dQ should be rather taken 
into account via the approximation of the domain f2 by domains with piecewise 
axes-aligned boundaries. Notice that for the anisotropic p— Laplace kind diffusions 

considered by Bendahmane, Karlsen in [17, LS], the discrete entropy inequalities 
on Cartesian meshes are as easy to obtain as for the isotropic case a(^) = k{£^)£^ 
considered in the present paper. 

• Taking into account sufficiently smooth dependencies on (t, x) of the convection 
and diffusion operators is possible, although quite technical; see [48, 11] for some 
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results in that direction, and also [33, GO] for well-posedness results for degener- 
ate equations with {t, x) dependent convection-diffusion operators. Discontinuous 
coefBcients are important for the modeling of fractured media. DDFV schemes 
for Leray-Lions operators div a(.T, Vw) with discontinuous (piecewise smooth) in x 
nonlinearity a are studied in the recent work [2:!]. The case of x-diseontinuous flux 
functions f{x, u) has received much attention in the last fifteen years (see, e.g., [llG] 
and the references cited therein), both from a theoretical and numerical perspective. 
Let us mention here that the problem of the choice of the appropriate entropy con- 
ditions strongly depends on the underlying physical interpretation; different models 
lead to qualitatively different admissible solutions. 

• Inhomogeneous Dirichlet boundary conditions can be taken into account, com- 
bining the techniques of [liT] with those of [1)]; both are rather involved, which 
explains our choice of the homogeneous boundary data for the presentation of the 
scheme and the convergence arguments. 

Appendix A: Proof of uniqueness 

This appendix is devoted to a proof of Theorem 2.2. The proof is an adaptation 
of the ones in Carrillo [-!)] (for entropy solutions) and that in Eymard, Gallouet, 
Herbin, Michel [45)] (for entropy process solutions, which can be viewed as entropy 
double-process solutions with /i = fj.*). The proof is mainly divided into several 
lemmas (Lemmas 9.1, 9.3, 9.5 below). For simplicity, let us only consider the case 
where the source term S is zero (see also Remark 9.1). 

We begin by introducing the set 

£ = {r e R : A^^{r) is neither empty nor a singleton}, 

and proving 

Lemma 9.1. Let {i',v*^v) be an entropy double-process solution of (1) with initial 
data VQ. Then for all W G R'', for any (j) e P([0, T)xn), c£R such that A{c) ^ £ 
and also for any (j) g I?([0,T) x 57), c G R^ such that A{c) ^ £, we have with the 
notation of Section 2 the following equality: 

i(7y±(^) + {d- l)rjt{^)) M + i( q±(z.) + (d - l)q±(z.*) ) • V0 
(76) - sign=^(t; - yl(c))(a(Vu) - W)) ■ V0 



dxdtda+ / ri^{vo)<p\t=o dx 



n 



eio 



lim / (signf )'(w - A(c))(a(Vt;) - VK) • Vu0 dxdt. 



Q 



Proof. We refer to [29, Lemma 1] and to [4!)] for details on the proof. The idea is to 
use tp :— (sign^(z; — A{c))(j)) as a test function in (D'.2). It is admissible; indeed, 
we can approximate it by functions in I?([0,T) x n) and pass to the limit in all 
terms of (D'.2), because V e L°°(g) nLP(0, T; W^'^{n)) for any of the two possible 
choices of (</>, c) (in particular, notice that sign^(w — A{c)) £ LP{0,T;Wq''^{^)) in 
case c G R''^). We have 

sign^(i/ — c) ~ sigii'^{v — A{c)) ~ sign^(z^* — c), 

thanks to the relation (DM) (which reads A{i') = v = A{i'*) in our notation) and 
to the choice of A{c) ^ £; then we use the weak chain rule to deal with the time 
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derivative. We also insert into (D'.2) the term 

' sign f{v- A{c))\/(f>-W - f f {signfYiv- A{c))Vw-W(l), 

JQ JQ JQ 

which is equal to = / div (M^sign^(w — A{c))(j)) for any of the two possible 

choices of (</>, c), by the Gauss-Green formula. As e | 0, the term containing 
(sign^)'(t; — A{c)) (f{v) — f(c)) • Vw vanishes, as shown in [29, Lemma 1]. □ 

We are now interested in comparing two entropy double-process solutions of 
(1), denoted by {v, ly* , 9) and {ji, fj* ^w)), of which the first one is chosen to satisfy 
V = V* . Consider the distribution I on 2?((5) defined by 



Jo JQ 



i((i/^/i)+ + (d-f)(^.-M*)+) dt^ 
+ i(sign+ (i.-/.) (fM-f(Ai)) 



(77) +(d-l)sign+(z.-A.*)(fH-f(At*))) ■ V0 

— sign''' {v~w) (a(Vw) — o(Vw)) • V(/) 

-I- / {vo — uo)^(l){0,x)dx. 



dx dt da dp 



Let us prove that we can write T as 
(78) I = IV+IM 

where TV[(j)] is defined by the analogue of (77) with each of v,v,vo, fi* ,w,uo 
replaced by its positive part; and is defined by the analogue of (77) with 

each of iy,v,vo, fx, ^* ,w,uo replaced by — — — Wq', — /i^, — (/i*)^, — w^, — Mq". 
To emphasize, whenever necessary, the dependency of J,JV,TAf on the involved 
solutions, we wiU write 1^]"^*^^, ,I'Pj;l;^'.%,IA/'J;l',^'"„, respectively. 

To justify (78), we use the identity (79) from the following easy lemma. 

Lemma 9.2. For all F : R ^ R such that F{0) ~ 0, for all a, 6 G R there holds 
sigTi + {a-b) (F(a)-F(6)) = sign + (a+ -6+) {F{a+)-F{b+)) 

+ sign + ((-a-)-(-6-)) iF{-a-)-Fi-b-)), 



(79) 



and 
(80) 



(i) sign- {b-a+) F{b) = -sign'^(a-^ - 6''') F{b+) + sign "(6) F(6), 
(m) sign"(fe- a-^) ^(0"^) = -sign '''(«+ - b+) F{a+). 



We apply (79) to a ~ b ~ ^ [or b ^ ^*) with F = Id and with F = f;, i = 
1, . . . , d. Futhermore, observe that the analogue of (79) still holds for a.e. {t, x) ^ Q 
if we take a = v{t,x), b = w(<, a;) and replace F{a),F{b) and F{±a^), F{±b^) 
by o(Vi;),a(Vw) and by o(±Vw''=),o(±V?i'''=), respectively. Indeed, we have, e.g., 
a( Vw) = a{Vv+) + a{\/v~) a.e. on Q, because v £ i''(0, T; Wq '^(17)). Using aU 
aforementioned identities, we split each term in the definition (77) of X into the 
sum of the corresponding terms in the definitions of XV and XJ\f. 

Now we estimate X "inside" the domain" . 
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Lemma 9.3. Let (i^, i/, w) and {^,,fj,*,w) be entropy double-process solutions of (1) 
with data vq, uq £ L°°{^), respectively. Then 1[4>\ > 0, G D([0, T) xQ), (/)> 0. 

Proof. The proof is an apphcation of the doubhng of variables method of Kruzhkov 
[03]; it follows [29, 49, 17]. We let ly depend on variables (t, a;,a) G Q x (0,1) 
and /i depend on another set of variables {s,y,/3) G Q x (0,1). In what follows, 
Vw means and Vw means VyW. As to the test function (/), it will depend 

on the variables {t,s,x,y), thus we will use the notations dt,ds and Va;,Vj, for 
the corresponding derivatives of (f>. We will work with nonnegative test functions 
4> £ T> (([0, T) X fi)^). Let us introduce the sets on which the diffusion term for the 
first, respectively, for the second solution degenerates: 

£^ = {{t, x)eQ I v{t, x) G £}, E^, = {(s, y) G Q I w(s, y) G £}. 

Denote by f^, the complementary sets in Q of fi/, f^, respectively. Observe that 
Vu = a.c. in and Vw = a.e. in E^ (recall (DM)). 

(i) First we apply Lemma 9.1 with the solution (y.,v,v). For all {s,y,i3) G 
E^ X (0, 1), choose W = a{\7w{s,y)) and take the entropy ?y^(-) = (• ~ c)+ with 
c = ^{s,y,f3), then with c = ^*{s,y,(3) in (76). We multiply the two resulting 
equations by ^ and by respectively, and add them together. Then we inte- 
grate in (s, y, P) e E^x (0, 1). Similarly, for (s, y, fi) & E^x (0, 1), we add together, 

with weights ^ and ^'^^"'"^ , respectively, the entropy inequalities (D'.3) for {i/,iy,v) 
corresponding to 77+(-) with c = fi{s,y,l3) and with c = i^i* {s , y , P) . We integrate 
the resulting inequality in {s,y,(3) & E^ x (0, 1). 

(ii) Next, we exchange the roles of {v^v,v) and {^,^*,w). This time we use 
the entropy ?7^(-) = (• — c)~; we use = a( Vu(i, x))\ and we only use one value 
c = J^(t,a;, a) in the analogue of (76) (for all {t,x,a) G E^ x (0,1)) and in the 
analogue of (D'.3) (for all {t,x,a) € E^ x (0,1)). 

(iii) Adding the inequalities obtained in (i),(ii), by the symmetry of the expres- 
sions involved (such as (y — /i)+ = {p, — i^)~ , etc.), we get, keeping in mind Remark 
2.1, the following inequality: 

fi /.I 



Jo 



QxC 



-{{u-^)+ + {d- 



l)[v~^i*)+){dt + ds 



sign+ im-Ki^)) 



+ {d- l)sign+ (f(i^)-f(Ai*)) )(V, + V, 



(81) 



■ sign+ (v-w) (a(Vw)-o(Vw)) • {V^ + Vy 
1 



dx dt dy ds da d(3 
{{vo-^i)+ + {d~ l){vo-^i*)+)(j)dx{dyds)dl3 



> lim 



(I' — uo)'^ (j){dx dt) dy da 
(sign+)'(v — w) (a(Vw) — a(Vw)) • (Vw — Vw) (j)dxdtdyds. 



The last term in (81) is nonnegative, because a is monotone and > 0. 

(iv) Let us now specify the test function. For /,n G N, let : R'' — ^ R, 
w' : R — > R be standard symmetric moUifiers with supports in {x G R'' | ||a;|| < M 
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and in {t G R I \t\ < j}, respectively. We take the test function in (81) to be 

0„j(i, X, S,y) ^(p {x, t) UJn {x - y) U)^ {t - s) = (j)Ulni^\ 

where G I?([0, T) x fi), > 0. With this choice, wc have 

(82) {dt + ds)<Pn,l = {dt(t>) (^nOj\ {V^ + Wy) (l)n,l = (V^,^) UJnOjK 

Then we let n,l — > oo. The first term in (81) converges to the first term in the 
right-hand side of (77). This argument is standard; one can use, e.g., the properties 
of the Lebesgue points of functions and the upper-semicontinuity of the "+- 
bracket" 



[uj], := / sign+(u) /+ / f 
Jn J{u=o} 



The two latter terms in the left-hand side of (81) are treated with the help of 
the triangular inequality and of the strong initial trace property (83) proved in 
Lemma 9.4 below. The limit, as n, ^ — > oo, of each of these terms is majoratcd by 
one half of the last term in (77) (this is because JgUj''{t)dt = ^ = /g ijj'(— s) ds). 
This concludes the proof of the lemma. □ 

Lemma 9.4. Let {fi,^*,w) be an entropy double-process solution of (1) with initial 
datum uq £ L°°(il). Then the initial datum is also taken in the following strong 
sense: 

Notice that another way to formulate (83) is to say that 

esslim / / ( ~,\l^^''^a\ -\ —\ii* — UQ\]dxda = 0, 

*io JnJo d J 

in the spirit of the original definition of Kruzhkov [63] . 

Proof. The proof follows the one of Panov in [7] , Proposition 1]. For c S R and 
/i > 0, consider the functions 

Ph{-;c): xefl^^J J (^^\fi(t,x]a)-c\ + ^^-^\fi*{t,x;a)~c\J dtda. 

Because ^* are bounded, the set {Ph{-',c)) is bounded in L°°(il). Therefore 
for any sequence /i„ 0, there exists a subsequence (not relabeled) such that for 
all c G Q, (pft (-; c) converges in L^{Q) weak-* to some limit denoted by p(-; c). 

Fix ^ G V{n), ^ > 0. From Definition 2.3, taking in (D'.2) test functions 
approaching ijj(t,x) := (l — ^)^^(x) we readily infer the inequalities 



Now for all S > 0, there exists a number N{S) G N, a collection {cf)fJ:P C 
and a partition of O into disjoint union of measurable sets fif , . . . , ^^^j-^) such th 
\\uo - UqIIli < (5, where 



(84) VcgQ / p{x; c) ^{x) dx < / \uo{x) ~ c\ ^{x) dx. 

Jn Jn 

By the density argument, we extend (84) to all ^ G L^{^), £, >0- 



N{5) 



Ely 
i=l 
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lim - — 



n Jo 



1, 



1, 



1/^ 



IM 



Un 



dt dx da 



= lim / > . ^ p,,„(a;;q)llf^5(a;)da; / V ^(2^; q )llf2* (2^) ^a; 



I (a;) da; 



Using once more the bound j|uo — WgHii < S (in the first term of the previous 
calculation), we can send S to zero and infer the analogue of (83), with a limit 
taken on some subsequence of {hn)n>i- Because {hn)n>i was an arbitrary sequence 
convergent to zero, (83) is justified. □ 

Lemma 9.3 tells us that X[-], which is a distribution on [0, T) x il, is nonnegative 
when restricted to 2?([0, T) x 57) and thus it is a locally finite measure on [0, T) x il. 

Now we show that X[(j}\ is nonnegative also for nonnegative test functions (f) that 
do not necessarily vanish on the boundary [0, T) x dfl. 

Lemma 9.5. Let [v^v^v), {fi,fi*^w) be entropy double-process solutions of (1) with 
initial data Vo,uq g L°°(yi), respectively. Then I[4i\ > 0, \f(f> G I?([0, T) x 57), > 0. 

Proof. We begin by modifying steps (i)-(iv) of the proof of the previous lemma; 
we refer to this proof for the notation and a part of the calculations. 

(i) We use (76) and (D'.3) in the same way as in the proof of Lemma 9.3; but we 
choose the values c = n~^{s,y,l3), c = {fj.*)~^{s,y, f3) and W = a(Vw+(s, j/)) instead 
of the values c = ^(s, y, /3), c ~ fi*{s, y, (3) and W = a(Vui(s, j/)), respectively. 

Notice that for all a, 6 S R, e > 0, we have sign j"(a — fe"*") = sign ^ (a"*" — 6+) and, 
moreover, this expression is zero whenever a < 0. Thus we can replace VjV, Vi; by 



,+ „+ 



Vw"^ everywhere in this calculation and obtain 



"'0 



QxQ 



+ (d- i)(z.+ -(/i*)+)+) 

+ i(sign+ (t.+ -/i+)(f(z.+ )-f(/i+)) 

+ (d- l)sign+ (z.+ -(M*)+) (f(^+)-f((M*)+)) )-V.^ 



(85) 



-sigii+ (o(Vt;+)-o(Vw+)) • V ^(t> 

1 



> 



dx dt dy ds da dfi 
(«-/i+)+ + {d-l){v+-{ii*)+)+)cpdx{dyds)dl3 
lim / / (signg!")'(w^ — {a{Vv^) — a{yw^))-Vv^(pdxdtdyds. 



(ii) We follow the proof of Lemma 9.3 but choose c = i'~^{t, x, a),W = a{\7v~^{t, x)) 
instead of c = v{t, x, a),W = a{\7w{t, x)). 

Let us apply identities (80) to a — v, b ~ 11 (or b — fi*) with = Id and with 
F = fi, i = 1, . . . , rf. Moreover, as in the proof of (78), we also have the analogue 
of (80) (i) with a = V, b = w with F{w) replaced by a( \7w). In the same way, we 
also have the analogue of (80)(m) with a ~ v, b ~ w, and F(v) replaced by a( Vu). 
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Furthermore, sign^(-) can be replaced by (sign^)'(-) in the above properties. In 
conclusion, we obtain 



Jo J JQxQ 



(86) 



+ -(sign+ (f(i^+)-f(Ai+)) 

+ (d- l)sign+ (f(;^+)-f((M*)+)) )-V,0 



■sign+ (t;+--u;+) (a(Vu+) -a(Vu;+)) • \7y(p 



dx dt dy ds da d(3 



U{v+ -u^y 4) {dxdt)dyda 
_ j X n 

> lim / / {sign+y {v+ ~w+) {a{Vw+)-a{\7v+)) ■Vw+ (f>dxdtdyds 



lim / / (sign^) {~w) aiV w) ■ V w (j) dx dt dy ds 



10 JO 



dgn-(/x){/ia,0+ (f(/i)-a(V«;)) • \7ycb} 



+ id~l)sign~{tl*){^i*^,(b+{Kp*)-a{yw))■ V, 
(mo)^ 4> ( da; dt) dy. 



dx dt dy ds da df3 



Qxn 



Notice that the sum of the last two terms in (86) can be rewritten under the form 
--^M.M'Cx), where 



(87) Xis,y)= / cf>{t,s,x,y)dtdxeVi[0,T)xn), 



and the distribution C^^^» is defined on 2?([0,T) x Q) by 



(x) 



{7Jo{^^)+id-l)vo{^^*))^sx 

JQ L" 



dydsd/3+ / r]„ {uo)xdy. 
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(iii) Adding (85) and (86), wc obtain, for any < e P (([0,r) x 17)^) with 
corresponding x defined in (87), the foUowing inequaUty: 



QxQ 



( + {d- !)(;.+ -(/i*)+)+) {dt + ds 



+ -(sign+ (f(^+)-f(A*+)) 

+ (d- l)sign+ K-(a.*)+) (f(^^+)-f((M*)+)) )-(V. + V, 



sign [V 



-) {a{\7v+)~a{\7w+)) ■ (V. + V^) 



dx dt dy ds da dfj 



(89) 



12 x( 



-M4 - + {d- m4 - {li'')+)+)c^dx{dyds) d(3 



JQxn 



[v'^ —Uq)'^ (p {dx dt) dy da 



> Um 
+ hm 

eiO 



(sign+)'(u+-w;+) (a(Vu+) -a(Vw+)) • (Vw+-Vu)+) c/ydxdtdyds 



(signj')'(— u;) a{\7w) ■ Vw (jjdxdtdyds — £^.^«(x) > —C. 



ix), 



where the last inequahty is due to the monotonicity of a(-). 

(iv) Now fix a;o G dil. Since dft is supposed sufficiently regular, there exists a 



vector Txf, and a positive number R^g such that the segment (cc, x - 



lies within 



il for all X G 90 n B{xo, Rxg), where B{x, R) stands for the ball of R"* with centre 
X and radius R. Choose in (89) the sequence of test functions 

2 



{t, x,s,y) = (l){y,s)ujn{ X ~ y + - 



LOl {t - s) = 4>U!niO 



for which (82) still holds. Notice that with this choice, the associated function 
Xn.i{y,s) in (87) writes as (l){y,s)9n{y)0\s), where 



On{y) ■■= / UJn{x-y + 



2 



dx. 



9'{s) 



uji{t ~ s) dt; 



n \\rxo\\. 

moreover, for all sufficiently large n e N we have 

cj^On e V{n) for aU e v([0,T) x (H n B(xo, i?.^))) ; 

9n{y) = 1 for all y e B{xa,Rxo) such that dist {y, 90) > |. 

As in the proof of Lemma 9.3, passing to the limit as Z,n — > cxd and taking into 
account the definition of the distribution TV, cf (78), from (89) we deduce 



(90) 



(91) 



^Vld] > -liminf 



/,n — >oo 

Now we remark that according to (D'.3), -C^,^* defined by (88) is a nonnegative 
distribution on [0, T) x O. Notice that the values of 6'„ are contained in the interval 
[0, 1]. Therefore for aU x G 25([0,T) xTl),(j)> 0, one has 

{<PSn] = C^,^, [0] - C^^^> [0(1 - 9^)] < £^,^. [</)]. 

It follows that 

(92) Z^,^. : X G V{[0, T) X H) ^ liminf (x0„) 
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is a nonnegative distribution on [0, T) x il; thus, it is a measure on [0, T) x fl. Since 
4> > and 6'' < 1, inequality (91) yields 

(93) IV[^] > - lim inf [xO^Oi] > - lim inf C^^^, [^9,,,] = -Z,,,^. [cf>] . 

Z.n^oo n— »oo 

It follows that XT' is a measure on [0, T) x fi. 

The remaining steps of the proof are aimed at showing, in an indirect way, that 
the positive part of the measure 2V does not charge the boundary [0, T) x dfl (in 
two particular cases, a direct proof of this fact is given in [73, 13]). Notice that 
this property is actually equivalent to the claim of the lemma; it accounts for the 
dissipative nature of the boundary condition imposed for entropy solutions. 

(v) Take (/> G 'd(^[Q,T) x {Ti D B{xo, R^„))y Fix m e N. It is easily checked 
from (90) that for all sufficiently large n G N, for all {t,x) G Q, 

0(s>y)(i - 0,n{y))On{y) = (l>{s,y)en{y) - (l>{s,y)em{y)- 

Therefore by (92), 

lim inf [(/)(1 - 0„i)] 

m— *C50 

= lim inf lim inf [0^n] — 1™ inf lim inf Cnn* [0^m] = 0- 
Applying (93) to the test fmiction 0(1 — O^n), we deduce 



(94) >limsupXK''''" 



(vi) Definition 2.3 of entropy double-process solution is invariant under the 
change of (/x, w),f,/,uo into (— /i, — /i*, — w),— f,— /,— uo. Moreover, one checks 
easily from the definition, cf. (78), that lA/'J^'.^'.^^ = ITZtlZf^^T . Therefore from 

(94) we deduce that for aU e d([0, T) x (H n B{xo,R^xo)f 

> lim sup X7'j;;^'"„[(/)6',„] +irZ'i-Zi^_l~'^[(l30,n] = hmsup2'^;^;"„[06'„] > 0, 

m — *oo L -I m — 'oo 

where the last inequality is due to (90) and Lemma 9.3. 

(vii) Not let (p be an arbitrary nonnegative function in 2?([0,T) x H). Choose 
a covering [j^^i B{xq, R^0, G N, of the compact set dfl. Introduce a partition 

of unity {£,i)iLo on ft associated with the covering U(UiI=i ^(^^O' ^xf,)) of ^^'^ 
apply Lemma 9.3 and the result of (vi) to the functions (/)^o £ 2-'([0,T) x il) and 
to e ^([O, T) X (Hn S(x^, i?^. ))) , i = 1, . . . , N, respectively. The claim of the 
lemma follows. □ 

Now we conclude the proof of Theorem 2.2. We have uq = vq. By a standard 
argument, choosing in Lemma 9.5 (j> = (f)(t) G I?([0,T)), we get for a.e. t G (0,T), 

(95) / / f^{iHt,x,a)-fi{t,x,P))+ + id-l){:,{t,x,a)-fi*{t,x,P))+)<0. 
Jo Jo Jn " 

Now, (95) means that for a.e. (x, a,/?) G f2 x (0, 1) x (0, 1), there holds 
ljL(t,x,P) = iy(t, x,a) — ^*(t,x,P), 



FV FOR DOUBLY NONLINEAR DEGENERATE EQUATIONS 



55 



which means that ^ — ^l* — v and each of them is independent of a, j3. This draws 
to a close the proof of Theorem 2.2. 

Remark 9.1. The proof of the contraction and comparison principle for entropy 
solutions of (1) (with S = 0) is essentially contained in the above proof. For nonzero 
source terms S, a more general version of inequalities (80) can be used; see [29] for 
the accurate treatment of this term. 



Appendix B: The reconstruction property 

Here we restate the result of [12, Lemma 8] and discuss its possible generaliza- 
tions. 

Lemma 9.6. Consider a triangle TT with vertices ^1,^2,^3 o.nd let to be the centre 
of its circumscribed circle. Denote by jlTj its area. For I G N/3N, denote by Ei 
the affine subspace < >; denote by TTi the triangle formed by toj 

and by {TTil its area, with the convention that the area is negative if to and ti lay on 
opposite sides from the line passing by Then 

2 ^ 

(96) ^El^'|P'^°jiJ.(^) = ^' foralir(^-K\ 

' ' 1=1 

Remark 9.2. For a multi-D generalization of the property (96), one could try to 
replace the projections on lines < > by projections on hyperplanes that contain 
the faces of the d-dimensional simplex TT. In this case one should replace the factor 
1^ by ^j3Y since |3r| = X]f=i^ l^'l ^^"^ because the dimension of Proj^j(T*) is 
[d — 1), whereas the dimension of is d. The proof of Lemma 9.6 given below 
shows that this generalization fails, except for very particular simplexes TT (this is 
clear from the multi-dimensional analogue of the identity (97) below). 

Remark 9.3. Using the "sine theorem", another proof of Lemma 9.6 can be given, 
which also works for any 2D polygon that admits a circumscribed circle. 

Proof. Proof of Lemma 9.6 For I G N/3N, denote by di the orthogonal projection if 
the point ti on the affine subspace Ei] set ~pi = todi and ~ai = toti. For /, i G N/3N, 
set b i^i = ~a i — ~ai. Denote by Tii the exterior to TT unit normal vector to Ei. 

Notice that we have for all / G N/3N, d i = [d i - it i)l% and also, for all i G N/3N 
such that i ^ I, 

_ ~pi ■ 111 

taking into account the sign of |2r;|. Since Proj^;^ +Proj^:^^^ is the identity opera- 
tor, (96) is equivalent to the statement that X]f=i l^dPi''3j< the identity 
operator. All vector V G R^ can be uniquely represented under the form 

3 3 
= ki ~a I with fc/ = 0, 
i=i 1=1 
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and thus, for all / G N/3N, V = Yl%i,i=i ( ~ "^0 = Yfi^i.i=i ^ Hence 

3 3 3 

A J2 |2r,|Proj<^,>(7^) = 2^ ^ 1?/ (T^ ■ = 2 ^ jtMh{b,, ■ iti) 

' ' 1=1 1 = 1 ' ' i,!^/; 1,1 = 1 ' ' 

33— 33 3 

We conclude that (96) is equivalent to the identity 

3 3 3 

(97) ki ~ai = —2 kij)i for all fci, . . . , /ca £ R such that fc/ = 0. 

i=i (=1 /=i 

Since to is the centre of the circumscribed circle of TT , the points di arc the centres 
of the corresponding segments Thus for all i,j € N/3N, by the Thales 

theorem we have p i — p j = —■^{ a i — a j). Hence (97) holds with ki G {0, 1, — 1}, 
j = 1, . . . , 3. Hence it holds for all choice of fc^. □ 

We refer to [7, 8] for a different kind of generalization of [12, Lemma 8] and a 
different proof of Lemma 9.6. 
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